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Abstract 

This is the theoretical review of exploration of new physics beyond the Standard Model 
(SM) in electric dipole moment (EDM) in elementary particles, atoms, and molecule. EDM 
is very important CP violating phenomenon and sensitive to new physics. 

Starting with the estimations of EDM of quarks-leptons in the SM, we explore the new 
signals beyond the SM. However, these works drive us to more wide fronteers where we 
serach fundamental physics using atoms and molecules and vice versa. 

Paramagnetic atoms and molecules have great enhancement factor on electron EDM. 
Diamagnetic atoms and molecules are very sensitive to nuclear P and T odd processes. 

Thus EDM becomes the key word not only of New Physics but also of unprecedented 
fruitful collaboration among particle, atomic, molecular physics. 

This review intends to help such collaboration over the wide range of physicists. 

I. INTRODUCTION 

This article is a review of the search of new physics beyond the Standard Model (SM) 
concentrating on electric dipole moments (EDM) of elementary particles like neutron, pro- 
ton, leptons, quarks as well as atoms and molecules. The presence of EDM implies T-odd 
and P-odd interactions. So if it exists, it indicates the direct T noninvariance as well as CP 
violation if CPT invariance is assumed. 

It is very important that these fundamental EDMs are enhanced in paramagnetic atoms 
{datom) and molecules (dmoiecuie) which have an unpaired electron. Also in diamagnetic atoms 
and molecules, proton and neutron EDMs appear via Schiff moment due to CP-violating 
hadron interactions. 

The discovery of CP violation in — )■ tt+tt^ decay in 1964 was an amazing event 
for the majority of theorists since the model at that time could not produce CP violation. 
The introduction of CP phase in the mixing matrix by Kobayashi-Maskawa was 7 years 
after that, which becomes the unique origin of CP violation in the SM jsl. This CP phase 
opened the door to new frontiers in a vast range of physics, especially in B factories: Belle 
at KEKB (KEK) and BaBar at PEP-II (SLAC). 
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CP violation in B mesons is measured by observing the asymmetry 
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Here B = B'^ = \db > or B = = \sb >, and / is a CP ei 



J/iJjKs, vr^vr , pKs- ^ means an imaginary part. For B'^ — ?• J/ipRg {4] |5|, 

S^A^/^ = 0.734 ±0.054, (3) 

with world average. It may be more advantageous in searching for new physics to consider 
the SM loop suppressed process like B^ — )■ (j)K^ etc. However, these results seem to be 



enstate such as 



consistent with the CKM mechanism {g]. 

As we will show, the EDM values predicted by the SM are very tiny because they appear 
first in three loops (quarks) and four loops (leptons) and are far smaller than the upper limit 
of the present and near future experiments. ^ On the other hand, there are some physical 
phenomena which suggest new physics beyond the SM other than neutrino oscillation ex- 
periments. The anomalous muon magnetic moment, = {g^ — 2)/2, is one such example 
8| 

- «r = (26-1 ± 8.0) X 10-1°, (4) 

corresponding to a 3.3cr discrepancy from the SM. There are also other indirect problems 
of the SM like the observed baryon asymmetry, n^/n^ = 1 x 10-^'^. Indeed, in the SM, CP 
violation is parametrized by the Jarlskog invariant, which is too tiny to produce this amount 
of asymmetry; we need other CP violating terms. The other implicit deficiencies of the SM 
are Dark Matter candidates and the hierarchy problem etc. 

Under these situations, the EDM is very important since some models beyond the SM 
give rather marginal predictions on the electron and neutron EDMs on the upper bounds of 
ongoing experiments. 



^ In the broad sense, there is another CP phase cahcd the 9 term in the SM, playing an essential role 
in especially the EDM of diamagnetic atoms (see section 5.3). 
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So new models are required to recover all such discrepancies. Furthermore, it must 
reproduce much larger phenomena which the SM predicts beautifully like, for instance, 
flavour changing neutral currents (FCNCs) and other vast low energy physics phenomena. 

Here we point out a peculiar property of the EDM: 

As is well known, EDMs of elementary particles are enhanced in atoms and molecules. 
In this sense, the EDM provides an unprecedented strategy of using atoms and molecules 
for the search of fundamental properties of elementary particles. 

Several review works on this subjects have been already published js] 10| 11| 12|. The 
new features of this review is that it is written by the author who is studying new physics 
beyond the SM, and, therefore, emphasis is on this point. However, EDM studies drive us 
necessarily to a wide range of physics (and chemistry), particle physics, atomic and molecular 
physics. The great achievements are possible only by the collaboration of theoretical and 
experimental scientists over this wide range of fields. Under these situations, we try in this 
review to make a small but significant bridge between these wide communities of scientists. 

Accordingly, we endeavor to give a self-complete concept of EDMs as far as possible, 
sometimes sacrificing the exhaustive citation of important references. 

II. BASICS OF EDM 

In this section we give the definitions and conventions used in this review, and basic 
formulae useful for the EDM. 

A. Definitions and Conventions 

Metric: 

.-.f' " V (5) 

\0 -l3x3 / 

Pauli matrices and spin matrices: 

-(;;)"■■(!:)"■(::) 

S' = (7) 
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Gamma matrices: 

l2x2 
— l2x2 

Chirality projection: 
Antisymmetric tensor: 
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-a* 



7^ = 27°7^7^7^ 




cr 



/it/ 
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-[7^7'^]^-(7V-7V). 



The electromagnetic field tensor is 



-e'^''B\ so = ^y^i 



IC. 



(8) 



(9) 



(10) 



The Cabbibo-Kobayashi-Maskawa matrix (2| and Jarlskog invariant |13| : 



V = 



C12C13 



'5l2Cl3 



Sue 



-iS 



"■512^23 ~ Ci2S23Si3e*^ C12C23 — Si2S23Si3e*^ S23C13 



S12S23 - Ci2C23Si3e 



i5 



-C12C23 - Si2S23Si3e C23C13 



(12) 



{u,C,i)YPLV 



d 

s 



(13) 
(14) 



JCP = \^iVajV^jV*,,Vi3k)\ = Si2S23Si3Ci2C23Ci3Sin(5. 

Here Jcp is the base independent CP phase called the Jarlskog parameter, appearing in 
CP violation processes via the Kobayashi-Maskawa mechanism. Hereafter, we denote the 
imaginary part (real part) of O by 5^(0) (JR:{0)). Apart from the EDM process discussed 
later, we also mention on neutrino oscillation processes. 



j<k 



(15) 



+ 4i ^ UajU*^jU*i,Upksm (ApjfcL) . 
j<k 

Thus we can determine the CP odd term (the third term) by measuring both P(z^/3 — ?■ I'a) 
and P{T'i3 — )■ T^n) - T2K [l^ found evidence of a nonzero 6*13 and recently the Daya Bay 
Collaboration 15| fixed it as 

sin22ei3 = 0.092 ± 0.016(stat) ± 0.005(syst). (16) 



Therefore the above mentioned CP phase experiments have become crucial. 



B. Effective Dipole Operator 



A permanent EDM of the electron must lie along its spin, namely d = dgcr [16 1. 
At tree level in the SM, a fermion ip (of mass and electric charge e, so electron's 
charge is e = — |e|) in the presence of electromagnetic field satisfies 

(7(p - eA) - m) ip = 0, ^ {'j{p + eA) + m) = 0. (17) 

However if we include loop corrections, the effective electromagnetic interaction is given in 
general by 

V = -eu2{p2)T^ui{pi) = J^A^ik) (18) 

with 

P = Pi+P2, k = p2-pi. (19) 

Here 

= (0, A) (20) 

is a true vector and transfoms as 

A^ (0, - A) under P, T transformation. (21) 

Whereas can be either a true or a pseudo vector. 

First we consider the case where is a true vector. In this case, takes the general 
form 

J'^ = Fi(A;2)(ii2ni)P^ + F2{k^)u2Yui + F3{e){u2U,)k>'. (22) 
However, from gauge invariance, the current is conserved 

k^J^' = (23) 

and 

Fsi^) = 0. (24) 
Using Gordon's decomposition for the bilinear form of a spinor of mass m, 

{u2a'"'ui) ky = -2mu2i^ui + U2UiP^, (25) 

the interaction term is then given by 



We should note that 

- = [S . B - Z75S • E] ^, 

where 

(27) 




Here off diagonal elements are suppressed by 0{v/c) relative to diagonal ones. The magnetic 
moment is defined as the coefficient of B in the above equation, 

= Q (^) cr, and /i„ = -2/1^, (28) 

where Q is the quark charge and m„ = md = 360 MeV has been assumed. It is clear from 
eq. fl27|) that the fermion has a magnetic dipole moment with g = 2 at tree level in the SM 
(!28|) . The second term is off-diagonal and there appears the additional P-odd term a'^p^ 
term in the product of the off diagonal element. The MDM and EDM of particles are defined 
in the rest frame and we hereafter neglect the off-diagonal element unless it is specified ^ 
On the other hand, a true axial vector current, the general form is 

In this case, survives due to chiral symmetry breaking. In weak interactions or higher 
loops in the SM or in new physics, the current includes both and J^'^ in general. Thus 
the following CP odd effective action appears, 

- V^^F^, = -^^ [zE ■ E - 755] ■ B] ^. 

4m 2m 

This will be discussed in more detail in connection with the EDM and MDM shortly. 
Also we can consider another conserved current like the vector case 

a{k^){'ykk^' - k^Yh^ (30) 

which reduces to 



^ This is true, especially for measuring EDM by spin precession as for most cases of neutral particles and 
atoms. It is not so serious for the measurements of EDMs of charged particles and neutral molecules 
though even in those cases, cooling down indicates a long coherent time is favourable. 
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in the nonrelativistic limit. 

This term is called aqn anapole term, which comes from the second term of 



.3 . -W) 



= / d'r'^^^, (32) 

in the expansion around r, that is, 

Af\v)= (vuVi-]T,ki (33) 



with 



T.M = -j dW,r[M'). (34) 

At the loop level in the SM and/or models beyond the SM, the following effective inter- 
action of gauge invariant form can be obtained: 



:^(A^- + 4)^a'^>F,. + l 

{Al + A)j)^E ■ + i{A% - Al)i}T, ■ Et/^.. (35) 



-{At + AWa^^^ljF,, + -{Al - Ai)^lja^^^'^F,^ 



Here we have neglected off diagonal parts in the second equalty. 

For the electric and magnetic dipole moment, we take zero momentum of the photon. 
Then the imaginary part of the coefficients of the effective interaction vanish because of the 
optical theorem (imaginary part of the forward scattering amplitude is given by the sum of 
possible cuts of intermediate states). We find the anomalous magnetic dipole moment 
and electric dipole moment to be 

= ^-^ = -^-^^{Al + Al), (36) 
= 2%Al-Al). (37) 

Note that Al and Ar must include a fermion mass (m^ or a fermion mass in the loop) 
because the effective interaction ipa^^ip changes the chirality which can be achieved by 
adding a mass term in the fundamental Lagrangian. If one of the particles in the loop is 
much heavier than the others, A^ and Ar are suppressed by the mass. Thus, for large A^ 
and/or Ar, it is preferred that masses of particles in the loop are similar to each other. 
The effective interaction ( 15^ also causes a £i — £^7 decay where the decay rate is given 

by 

777 

m-^^n) = ^{\^t? + \^l?) ■ (38) 
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Thus EDM and MDM have opposite parities and different on the order of magnitude. How- 
ever they appears in parallel, and have some similarities also. One of them is their SU(6) 



property [17| and will be discussed in Appendix A. 

For an invariant electromagnetic field, EDM, MDM, anapole, and higher n-pole moments 
appear as the multipole expansions of the Coulomb potential and vector potential. These 
points are also discussed in Appendix B. 

Quarks receive additional contributions, which will be discussed for diamagnetic atoms. 
Here we list the results. 

A strong CP violating term connected with the 9 vacuum (see Appendix G) 

L,=4 = -^'OGl^Gl^e^'^^^ ^ :^dG^ ■ G\ (39) 

In new physics beyond the SM, we have the other P and T violating effective actions: 
the chromoelectric dipole operator (cEDM) 

Lc = -'-d,qa,,^,TyG^'"^ = -'-d.qaGj^q, (40) 

and the following dimension 6 operators, 

Lg = ~-^dGfabcG'^^pG''^^G1^e^''^'^ = —-dcfabcG'^G^G^, (41) 



the so-called Weinberg term 



18|, and 



Ld=6 = J2 C^,^P^Oa^piiJ,Oal,i^r (42) 



Here ipi and ipj are leptons and/or nucleons. Oa are scalar, vector, and tensor gamma 
matrices. We will explain the detailed physical implications in the diamagnetic atom in 
section 6. 

In the SM, weak interactions act with matter and gauge fields in the form, 

H^eak = i^PiT^W = Jf.W''- (43) 

However, except for the top quark, fermion masses are small compared to weak bosons 
masses. They are descrived as the four Fermions coupling 

Hyjeak J/iJ^- (44) 

This is the case for tree diagrams. If you consider loop diagrams or new physics beyond the 
SM, we will encounter more general forms. We will discuss this in the Appendix C. 



C. Experimental Bounds 



We have no experimental signal of the EDM yet but have upper limits. They are as 



follows 



19|. It is very impressive that recently we have a more precise upper limit of the 



electron EDM from molecule (YBF) than from atom (Tl). 



4 from thallium atom d{Tl) = (6.9 ± 7.4) x 10"^** e cm j20| 



df, = (3.7 ±3.4) X 10-^^ecm (21| 



dn < 2.9 X 10"^^ e cm (90%C.L. 
dC^^Hg) < 3.1 X 10-29 e cm {95%C.L. 
de from the molecule d{Yh¥) = (-2.4 ± 5.7stat ± 1-5 



22| 



23| 



'syst ) 



X 10-2^e cm 



24| 



(45) 
(46) 
(47) 
(48) 
(49) 



For reference, we give here the muon anomalous MDM, a non-null signal of new physics 
beyond the SM. 

The deviations of the SM predictions from the experimental result are given by 



cxp 
cxp 



14.8(8.2) X 10"^°, 
] = 30.3(8.1) X 10"^°. 



(50) 



Here the hadronic contributions to [t] and [e"'"e ] were calculated [25| by using data of 



hadronic r decay and e^e~ annihilation to hadrons, respectively. These values of Aa^[r] and 
Aa^j[e~'"e"] correspond to 1.8(T and 3.7(T deviations from SM predictions, respectively. EDMs 
and MDMs comes from similar diagrams apart from CP transfomation and the differences of 
magnitudes in the MDM and EDM stems from the cancellation of diagrams and symmetry. 
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III. STANDARD MODEL 



In this section we give the EDMs of quarks, hadrons, and leptons in the SM framework. 
Structure of matter multiplets in SM+(Dirac) neutrino is 

= (ul u'i ul) ~ (^3, 1, ^ 

3 



d'==K dl dl)r^ Ul,^), (51) 




CP violation occurs a la Kobayashi-Maskawa mechanism, that is, CP phase in CKM 
mixing matrix for quarks or MNS matrix for leptons. Diagrammatically, it resembles with 
LFV processes but the former is necessary to incorporate the non-zero Jarlskog parameter. 
Apart from the uncovered new phenomena in neutrino, the SM has the deficiency of baryon 



assymmetry. Jarlskog introduced Acp 



131 defined by 



[M^MlMdMl]=iAcp. (52) 

Here M^, are up- type and down- type quark mass matrices. The observables are not 
these matrices but those which are invariant under rebasing and rephasing, that is, eigen 
values and CKM mixing matrix. By construction, Aqp is traceless and Hermitian, and 
characterizes the effect of CP violation. Its explicit form is 

det Acp = [m] - ml){ml - ml){ml - ml){ml - ml){ml - m^)(m^ - ml)Jcp, (53) 

where Jcp is given by f|T^ . As we will see the detail shortly, one-loop diagram (Fig. ([1]) 
gives zero contribution to the EDM. As for second loop (Fig. |2]), using this J and denoting 



by / the Green function of / flavoured fermion 



261], the f quark EDM has the form 



J2 %v,,VifV*jV;,)fjkif = I'^iv^kVuV/fViDfUki - ikj). (54) 



jkl 
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FIG. 1: The diagram for the EDM of d quark at one loop level in the SM. 

and J ^ 3 X 10~^is twice the area of unitary triangle V*^Vub + Vcd^cb + VtWtd- One finds 
that J/T^y^ < 10^^'^. This falls short of the baryon asymmetry ns/n^ by 10""'^°. 
We will show the detail loop by loop in the subsequent sections. 

A. Quark EDM 

For definiteness, we consider the EDM of d quark. The advantage of d in comparison with 
u may be that t can be in the loop with the weak interaction to avoid the GIM cancellation 



27| 



1. One Loop 

In the one loop level (Fig. [T]), elements of the CKM matrix appear as iVg^p. Therefore, 
there is no EDM (imaginary part of coefficient) apparently. 

2. Two Loop 

At two loop level in the SM, two types of diagrams have an imaginary coefficient poten- 
tially. The diagrams are shown in Fig. [2l It is, however, shown that the imaginary part 
of each diagram vanishes by the summation of contributions from all quarks of internal 
lines 



For the diagrams in Fig. [21 it is clear that q" must not be d quark because it gives 
l^g'dPl^'gdp of real value. By the same reason, q' must be different from q. 
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FIG. 3: Diagrams for the EDM of d quark at three loop level in the SM 29|]. 



3. Three Loops 



Fig. [3] shows three loops diagrams which contribute to d quark EDM in the SM. The 



formula for the contribution is given in 



29| as 



e ~ IOSttS 



^bc ~ '^^bc + ] Lwb + —Lie 



335 2 



1231 7 



36 ■rl'--l0S^r^'<^'^ 



(55) 



where Lab = lii('^a/'^6)- results in 



dri ^ —10 e cm 



(56) 



while the triple log approximation (taking only term) gives 



dd ~ +10 ^^e cm. 



(57) 
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FIG. 4: Diagrams for the EDM of neutron in the SM 



3C|. 



B. Neutron EDM 

The dominant contribution to the neutron EDM in the SM comes from "two loop" dia- 
gram in Fig. m The interaction on the left part of the loop is given by the phenomenological 
interaction hamiltonian 

H = iGFmlun{A + B-f^)uj:^^, (58) 
A = -1.93, B = -0.65, (59) 

where m„, Ms, and (/?^ stand for wave functions of neutron, E~ baryon (dds), and tt"*", 
respectively. The interaction on the right part of the loop is so-called "strong penguin" 
whose effective operator is given by 

12v2 7r mi 

^ ^ q=u,d 

Note that the "Hpcn seems to be obtained for rrtc c:^ rrit < my/. With these interactions, the 
neutron EDM was estimated as 

= df°'' + ^ 10"32 e cm . (61) 

Here the first is the contribution from Fig. [3] (©(agGp) ~ 10~^^ ecm) and the second from 



Fig. 1 130| 



If we incorporate the rephasing invariance of strange wave function, this value is modified 
to 1.4 X 10^3^ < \dn\ < 9.9 X 10-33 e cm |31|. 



There are new CP violating five and six dimensional operators, 

Ccpv = d,q{cTF)^^q + d,q{cTGh,q + wGGG + ... (62) 
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FIG. 5: the EDM as a function of the pseudoscalar meson mass squared m'j:,^ for neutrons (left 
panel) and protons (right panel). The arrow shows the physical point of the pion mass squared, 



0.0195GeV^, and the star symbol denotes the the result of the current algebra (C.A.) 



32|. 
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W W 
FIG. 6: Diagram for the EDM of W boson at two loop in the SM. 



The detail of this contribution will be discussed in the section of diamagnetic atom. 

The concrete and model independent calculations are expected in the lattice QCD. We 
will list only their results. It is also possible to guess the rough estimate of hadron EDM 
using SU(6). The detailed discussion is given in Appendix A. The EDM of neutron is 
estimated by lattice calculation (see Fig. [5]). 

C. Lepton EDM 

For definitteness, let us concentrate on electron. Similarly to quark EDM, one and two 
loop diagrams do not contribute to electron EDM. In order to avoid GIM cancelation 
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W W 
FIG. 7: Three loop diagram which may give a nonzero contribution to the EDM of W boson. 



(oc (m^ — m|)/m^), CKM matrix is better to be used than Maki-Nakagawa-Sakata matrix 
(lepton mixing). Two W bosons (at least) should attached to the electron line in order to 
use a quark loop. Then, the electron EDM is caused by the W boson EDM. It was shown 



that the W boson EDM vanishes at two loop level 33|. 
The two loop diagram is shown in Fig. [6l 

Jcp defined by (fT4|) is antisymmetric under j f-)- / (IMI) (corresponding to side line's 
quarks), whereas it is symmetric in Fig. 5. Adding aother loop with gluon (see Fig. [7]), the 
W boson EDM in three loop was estimated as 

The electron EDM in four loop was estimated with dw as 

2 

4 ^ J—^H^ dw^Sx 10"^^ e cm . (64) 
32%^ mw 

IV. BEYOND THE STANDARD MODEL 



In this section we discuss models beyond the SM. 

As we have shown in the previous section, the SM predicts rather smaller values of EDMs 
than the experimental upper limits by roughly ten orders of magnitude. However, we also 
know that CP violation in the SM is insufficient for baryon asymmetry in the real world. 
Also we have many direct signals of new physics beyond the SM like neutrino oscillations 
and muon g-2 etc. Even if we stand in the SM, we have new 4- and 6-dimensional CP 
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violating effective actions like (139|) to (H2|) . which have never been discussed so much in the 
previous section. Then it is very natural to estimate how much such new physics or new 
models predict the EDMs. In this section we concentrate on new physics beyond the SM. 
As for the new 4- and 6-dimensional CP violating effective actions, we will discuss in the 
sections of diamagnetic atoms and molecules. 



A. Minimal Supersymmetric Standard Model 

In the MSSM, all particles have their SUSY partners; sfermions / (bosons) for fermions 
/s, Higgsinos (fermions) for Higgs bosons, gauginos (fermions) for gauge bosons, and and 
another Higgs doublet is added to the SM for recovering the chiral anomaly free condition 
once broken by this doubling js^. Yukawa coupling is given by 

Wmssm = YuuQH^ - Y^dQHa - Y^eLHa + ^iH^H^ (65) 

SUSY is broken at O(lTeV) by the soft SUSY breaking term which retain hierarchy problem. 
MSSM is the minimally extended supersymmetric SM and we will consider below cMSSM 
and i^MSSM including light neutrino mass. In general soft breaking terms are 

+ f^lA^eR, + filHlH, + f,lHlH2 + filS*S (66) 
+ ^uijU^aQuHi + Adijd*^iQLiH2 + Aeije}^iQLiH2 

ifj^ = 1 i^M^ig^'g" + M^W^'W'' + MiB^'B" + c.c.) (67) 

These include many CP violating phases, in general. 

For the loop correction to the Higgs masses, the problem on its quadratic divergence is 
cancelled by the loop of SUSY partners (different statistics with same coupling). 

SUSY particles contribute to fermion EDM at one loop shown in Fig. [HI The neutralinos 
X° and charginos are mass eigenstates, and they are linear combinations of Higgsinos 



and gauginos of SU{2)l and U{1)y- The d quark EDM from the diagram is estimated |35| 
as 

, / 1 . Ix4«. ^ , f^m (I ^ ^ /i^ /i^(/x^ + 2M^) , M2\ 

(68) 
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Ir fh 

/ \ 
/ \ 

Ir < I < I < fl 

FIG. 8: Diagrams of SUSY contributions to the EDM of fermions. Gluino g contributes only for 
quark EDM. 




Here v is the common vacuum expectation value of Hi and H2. m and M are masses 
of d quark and the universal squark mass (by assumtion), respectively. If we adopt 
v'^iV^^AoV^) ~ M2 with rough estimations of M = 100 GeV and maximal mixings, we 
obtain 

10-22 e cm. (69) 

The value is clearly in conflict with the experimental bound on dn- The naive estimation 
was done with Mq ~ 100 GeV and a sizable CP- violating phase sin0 ~ 1. Therefore, the 
contradiction can be resolved by small phase (approximate CP symmetry) and/or heavy 
masses of SUSY particles. However, we adopt not such fine tuning but the universal soft 
SUSY breaking (constrained MSSM; cMSSM). Vl{r) is the unitary matrix which rotate left 
(right)-handed weak eigen state, and therefore if Aij oc Yij, V^^A^Vj^ is real diagonal and 
the imaginary part of it vanishes. Thus the small EDM leads us to the relations of trilinear 
terms in (165|) and fl67|) 

A^ = A^oyu, Aa = AdoYd, A, = A,qY,. (70) 

Also we accept the universal soft SUSY breaking which is realized by gravity or gauge 
mediated SUSY breaking. 

Mi,M2,M3 ~mi/2 (71) 
mg, rriL, niu, rrid, nie, mn^ , m^a ~ ruQ. (72) 

CP-violating phases appear only in flavor off-diagonal parts of matrices (hermitian Yukawa 
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FIG. 9: Diagrams of SUSY contributions to the EDM of fermions in case of flavored CP violation 
where CP violating phase appears with change of flavor 36 1. 



matrices) and the CP violating effect is suppressed by small mixings only due to RGE. 



(m 



2 ' "RR 



[m. 



m 



, 5 



RR 



m 



2 ' 



m 



2 ' 



JRR 



(73) 



The diagrams for the flavored case are shown in Fig. [91 When both left- and right-handed 
squarks (sleptons) have mixings, they contribute to the EDM in the form: 

Jlr = ^Ay^iSlLh, J^R = ^{S'-RRVuSlLh. (74) 



It was shown (see e.g. 36j) that d^^ can be ~ 10"^^-10~^^e cm, and then SUSY parameters 
are constrained by hadronic EDM. 

Also there are additional diagrams called Barr-Zee diagrams that contribute to the EDM 
beyond one loop level (see Fig JTOi) . 

Thus the MSSM gives an elegant backgroud but itself does not predict any definite 
relations between quarks and leptons including all observations in neutrino. 

In other word, its predictions are not testified from many constraints from various already 
known observations. These observations must be complicatedly related in reliable models. 

It is a Grand Unified Theory (GUT) which fulfiles these deficiencies 39 1. 



B. Minimal Supersymmetric SU(5) GUT 



Here and hereafter "minimal" means the minimum number of Higgs field with renormal- 
izable Yukawa coupling. 
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7, (S) 




FIG. 10: Barr-Zee diagrams: the Hi lines denote all three neutral Higgs bosons, including CP- 
violating Higgs-boson mixing, and heavy dots indicate resummation of threshold corrections to the 
corresponding Yukawa couplings js?]. 



In SU(5) model [38|, matter multiplets in ( 152|) are classified to 



5* 



e 

V ^ / 



10 



/ — ^2 Ml di ^ 

Ul U2 d2 

Us d-i 
e+ 

\ J 



(75) 



We need two Higgs h*^ and 24^, and SU(5) breaks down to SU{'i)c x SU{2)l x f/(l)y by 

24.H = diag {v, V, V, -h^, . (76) 

Here 5 and 10 are broken to 



5 = (1, 2)(l/2) + (3, l)(-l/3)), 10 = (1, 1)(1) + (3, l)(-2/3) + (3, 2)(l/6). 
Then, SU{?>)c x SU{2)l x f/(l)y breaks down to 5f/(3) x U{1)q via 

5* = (0,0,0,0,t;/y2). 

Yukawa coupling has the form, 
W 



^/^.10a0,5^ + v^/^.10,5*,5J, + /^-S^I.-Sh + M,,ia 



r 



(77) 



(78) 



(79) 
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Here the products imply 

10,10,5^ = eabcdelOflOfblj 

10i5*j5*jj = lOf 5f 5^* etc. (80) 
with a,b = 1, 5. Then mass matrices have the following forms 

M'^ = = fv/V2, = rv/V2 (81) 



at GUT scale. This gives nice b — t unification. The disparity between their observed 
masses is supposed to be due to renormalization effect from Mqut fo their mass shell. 
Unfortunately, we can not explain the disparities between the first and second families even 
if we renormalization effects since it predicts wrong relation 

md/ms = me/m^, mjmi, = m^/m^. (82) 

t also predicts too fast proton decay 



4l|. 



401 1 . Hadronic EDM in SUSY SU(5) was discussed in 



Flipped SU(5) changes 

u^^d", e^^u" (83) 
and, therefore, we obtain in place of ( 18T|) 

M„ = M,. (84) 

This does not leads to apparent pathology. Moreover, it is attractive from Doublet-Triplet 
splitting: Higgs super potential has the form 

WO? = 10 X 10 X 5 + To X To X 5, (85) 

which give rise to triplet mass 

((1, 1, 0)io)(3, 1; l/3)io(3, 1 : -1/3)5 + ((1, 1; 0)to)(3, 1; -l/3)io(3, 1; 1/3)5 (§6) 

but has no doublet mass since 5 + 5 has no partner in 10 + 10 (Missing partner mechanism). 
This is a solution to the Doulet-Triplet problem without additional adjoint Higgs. However, 
flipped SU(5) drives us to unrenormalizable heavy Majorana neutrino mass term, 

lOilOjTOj^TOj^/A (87) 
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for seesaw mechanism. 

The other approaches are to introduce unrenormahzable term 



Wy = e,fe,,e(/wl0fl0f24|,_^5{, + /2zjl0fl0f24^5^)/A 

+ guj5*H,24aHa0T5*JA + to,5l,,10f 24^,5*,/A (88) 

or to add another Higgs in Yukawa couphng, 

Y,,5*10,45% (89) 

etc. Unfortunately in SU(5) model, right-handed heavy Majorana neutrino belongs to the 
singlet and we have no constraint on it. Usually it is assumed to be diagonal but there 
is no reason to be justified for that. The other undetermined parameter dependence (like 
niQ. Ml/2, Aq, tan/3) crucially depends on this assumption. 

These points are remedied in case of renormalizable SO(IO) GUT, which is discussed in 
the next subsection (cEDM and parity odd nuclear interaction). 

C. Minimal Supersymmetric SO(IO) GUT 



In the SO(IO) Grand Unified Theory 43|], fermions belong to a multiplet of 16 represen- 
tation as 

ijj = (m^, wfj, M^, rf^, d% rf^, cr, ur, ul, ul, ul, dl, dl, d^, ei, t^l,)^ ■ (90) 
Note that the right-handed neutrino ur is included naturally. 



So-called minimal renormalizable SO(IO) model includes Higgs bosons of 10 and 126 in 
Yukawa couplings. This is because 

16 x 16 = 10 + 120 + 126. (91) 

In order to make singlet in Yukawa renormalizable coupling, therefore, Higgs can be 



10, 120, 126. ^ One Yukawa coupling leads to the conclusion that the CKM mass ma- 



trix is unity and we needs at least (minimal) two Higgs, 10 + 120 or 10 + 126. We select 



^ If we relax the renormalizability , different SO(IO) models are also possible 4^. However, in this case, 
we have much less prcdictivity 
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the latter set. This is because 



126 = (6, 1, 1) + (10, 1,3) + (10, 3, 1) + (15, 2, 2) (92) 

under SU (4)^ x SU (2) ^ x SU {2)ji and the second and third terms play essential role in type 

45| . 126 is necessary to be added. 



I and type II seesaw, respectively. In its SUSY version 
Providing the Higgs VEVs, Hu = v sin /3 and Ha = v cos /3 with v = 174GeV, the quark and 
lepton mass matrices can be read off as 

Mu = CioMio + Ci26Mi26 

Md = Mio + Mi26 

Md = CloMio - 3ci26Mi26 (93) 

Me = Mio - 3Mi26 
Mt = CtMi2q 
Mr = CrMi26 , 

where M„, M^, Md, Me, Mt, and Mr denote the up-type quarlc, down-type quark, Dirac 
neutrino, charged-lepton, left-handed Majorana, and right-handed Majorana neutrino mass 
matrices, respectively. Note that all the quark and lepton mass matrices are characterized 
by only two basic mass matrices, Mio and M126, and four complex coefficients cio, Ci26, ct 
and cr, which are defined as Mio = Yi^a'^'v cos (3 , Mue = Yue f3'^v cos f3, cio = (a"/a;'') tan/3, 
Ci26 = (Z?"//?*^) tan/3, Ct = vt / {P'^v cos (3)) and cr = vr/ {f3'^v cos (3)), respectively. These 
are the mass matrix relations required by the minimal SO(IO) model. The model is very 
predictive by virtue of the relation between quark Yukawa matrix, lepton Yukawa matrix. 



and neutrino Majorana matrix. Fig. [TT] 45| shows the prediction for the electron EDM \de\ 
in the minimal SUSY SO(IO) with respect to the universal gaugino mass Mi/2- The muon 
EDM \df^\ exists above \de\ by roughly a factor of 10^. The muon anomalous MDM and 
the decay branching ratio of /i — )■ 67 (Lepton Flavour Violation: LFV) are also predicted 
(see Fig. [TTl) . EDM anf LFV are due to the essentially same diagrams, apart from the fact 
that the former (latter) comes from diagonal (offdiagonal) part of sfermion mass matrix.^ 



^ We have added the last panel with nonzero Aq reflecting the recent discovery of Higgs-like particle around 
125 GeV by the LHC. See the last part of Discussion for more detail. 
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FIG. 11: The predictions for the electron EDM \de\, the muon anomalous MDM (jSOp . and 
the decay branching ratio of — )• 67 in the minimal SUSY SO(IO) with respect to the universal 
gaugino mass Mx/2- Trilinear term ^0 is assumed to be zero except for the last panel. The last 
panel is added for the reference to see the behaviour of non zero Aq, where the branching ratios, 
Br(T —7- /i7) (top) and Br(^ — t- 67) (bottom) are given as functions of Aq (GeV) for mg = 600 GeV 



and Ml/2 = 800 GeV. Ah are cited from 



i5|. 



The effective Lagrangian relevant for the EDM, MDM, and the LFV processes {ii — ?• ij'y) 
is described in (135|) . R, L = {1 ± 75)/2 is the chirality projection 

.e 



^eff = -^^rne-i,a,:.F^- {A'^Pl + A^^Pr) 



(94) 



where Pl,r are Left-Right projection operators, and Al^r the photon-penguin couplings of 1- 
loop diagrams in which chargino-sneutrino and neutralino-charged slepton are propagating. 



by ^. 



The 



it should be noted that we have changed the normalization of Ai^^ji from 
explicit formulas of Al^r etc. used in our analysis are summarized in 45|] j46|. If the diagonal 
components of Al^ji have imaginary parts, the EDMs of the charged leptons are given by 



dije = -mi^Q{A 



A 



R) 



in the new normalization. The rate of the LFV decay of charged-leptons is given by 



IGtt 



-rrio 



{\A'l? + \A=R?) 



(95) 



(96) 
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while the real diagonal components of Al^r contribute to the anomalous magnetic moments 
of the charged-leptons such as 



6a 



SUSY 



-m 



(97) 



In order to clarify the parameter dependence of the decay amplitude, we give here an ap- 
proximate formula of the LFV decay rate 46|, 

' ^ |2 



r(£. ^ £,7) 



-m», X 



a2 



Ml 



tan^ 13 



(98) 



where Ms is the average slepton mass at the electroweak scale, and ^ Am| j is the slepton 
mass estimated in Eq. ( l99l) . We can see that the neutrino Dirac Yukawa coupling matrix 
plays the crucial role in calculations of the LFV processes. We use the neutrino Dirac Yukawa 
coupling matrix of Eq. fllOOp in our numerical calculations. In the leading-logarithmic ap- 
proximation, the off-diagonal components (i 7^ j) of the left-handed slepton mass matrix are 
estimated as 



3m2 + Al 



(99) 



where the distinct thresholds of the right-handed Majorana neutrinos are taken into account 
by the matrix L = log[MG / Mfi-]6ij . 

Unlike the muon MDM, quark and lepton EDM has still null observation. This is of 
course due to tiny CP violation and due to the cancellation of diagrams where 7 (gluon) 
couple with slepton (squark) and where they do with Higgsino (gluino) in (Figj9]) 47 |. 

If we consider gauge mediation scenario for SUSY breaking, Aq ^ 0. In the basis where 
both of the charged- lepton and right-handed Majorana neutrino mass matrices are diagonal 
with real and positive eigenvalues, the neutrino Dirac Yukawa coupling matrix at the GUT 
scale is found to be 
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-0.000135 - 0.00273Z 0.00113 + 0.0136z 0.0339 + 0.0580z ^ 



v 



0.00759 + 0.0119i -0.0270 - 0.00419i -0.272 - 0.175i 
-0.0280 + 0.00397^ 0.0635 - 0.0119i 0.491 - 0.526^ 



(100) 



Semi-Leptonic LFV processes was discussed in 49 1. 

Thus the EDM, MDM, lepton flavor violations etc. are all closely connected, which are ex- 
pected to be explained universally by GUT. However, we do not adhere to the special model 
in this review and also discuss more phenomenological models in the following subsections. 
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These may be the remnants from GUT or may be independent on GUT. For instance, 
the adjoint representation of SO(IO), 45 is decomposed into 

45 = (1, 3, 1) + (3, 1, 1) + (15, 1, 1) + (6, 2, 2) (101) 

under S'f/(4)c x SU{2)l x SU{2)r and leads to Left-Right symmetric model, ql = Qr- Also 
10 representation is decomposed into 

10 =(1,2, 2) + (6, 1,1), (102) 

which leads us to two Higgs SU{2)l doublets under the SM. Also 126 is 

T26 = (6, 1, 1) + (10, 3, 1) + (To, 1, 3) + (15, 2, 2). (103) 

If (10, 1, 3) ((10, 3, 1)) has vev, it gives type I (type II, or Higgs triplet Model) seesaw model. 
In the following we consider these models independently on GUT principally. ^ 



D. Two Higgs Doublet Model 

Most of models beyond the SM has some new Higgs bosons. As the simplest extension of 
the Higgs sector of the SM which has only one Higgs doublet Hi, another Higgs doublet H2 



is introduced in the Two Higgs Doublet Model 



5l| . There are several types of the model 



depending on which doublet couples with which fermion: 

type I (SM-like) : Hi couples with all fermions 
H2 decouples with fermions 
type II (MSSM-like) : Hi couples with down-type quarks and charged leptons 

H2 couples with up-type quarks 
type HI (general) : both of Higgs doublets couple with all fermions 
etc. etc. 

If CP-violating term exists in the Higgs potential, e.g. [H\Hi)[H\H2) with an imaginary 
coefficient, there appears the mixing between CP-even {H^) and CP-odd {A^) neutral Higgs 



^ He et al. discussed neutron EDM in those models 



50|. 
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FIG. 12: Diagrams of muon EDM in two Higgs doublet model [52i] . 



bosons. Then these Higgs bosons can contribute to the EDM (see Fig JT2|) . The mix- 
ing between and A° provides also CP-violating electron-nucleon effective interactions 
[ei'-f^eNN , etc.) which will contribute to the atomic EDM. 



Barger et al. 



52| gave large close to the propose experiments. It should be remarked 



that their values are calculated in units of '^Z where 



Amplitude(g^) = 5^-^^?^, (104) 



n 



and it is probable that \'^Z\ ^ 0. Indeed, the masses of neutral and charged Higgses 
and phases are tightly constrained from Rb = r{z^h<ldrons) ' ~^ g7), — B mixing, p 
parameter etc., and we should take those constraints into account 



53| 



E. Higgs Triplet Model 



In the Higgs Triplet Model |5J], we introduce a SU(2) triplet Y = 2 scalar as 
A+/V2 A++ \ _ 2 

^ J ^tripletYukawa = —hasL^ia ^PlLr + h.C. (105) 

AO -A+/V2 J 

This model generates neutrino masses without right-handed neutrinos with the triplet vac- 
uum expectation value V/\ which is given by the explicit breaking of the lepton number. 
This model is very predictive because of a clear relation 

niap = V2VAhal3, (106) 

27 



where m^/? denotes the Majorana mass matrix for neutrinos. There is no new interaction 
with quarks and no effect on quark EDM. Unfortunately, this model can not give a large 
contribution to lepton EDM also because of the absence of the new interaction with right- 
handed fermions. For example, one loop diagram for electron has a factor of |/iaeP (similarly 
to Fig. [I]). 



F. Left-Right Symmetric Model 



Left-Right model 55j is used in a varieties of ways and needed to be clarified. If we 
consider it as a remnants from SO(IO), 5*0(10) SU{A)c x SU(2)l x SU(2)r, it satisfies 
at vps energy scale 

9l = 9r (107) 

and PS model is unified at Mqut as 

M4 M2L M2LR Ml/2 

= = = . [iUii) 

«4 0!2L a2R OiGUT 

Also mixing matrices of left-handed and right-handed fermions are same. Of course these 
constraints are realized at vps but deviate from as energy goes down to the SM scale by 
renormalization effects. 

However, if we consider a model of SU (2) l x SU (2) nxU (1)b_L) we are free from the above 
constraints. For instance, in the framework of SO(IO) GUT, vr is order of O(10^^)GeV. 
However, if we go apart from GUT but still consider that the mixing matrix of right-handed 
quarks Vr has a similar structure as that of left-handed quarks Vl, the lower limit of 
is relaxed to Mwr > 1-6 TeV 56|. Moreover we may go beyond restricting on V-A and 



V-l-A interactions only and may consider general form; 



9RR{eR^eL)0^Lf^R) + gRL{eR^eL)i'i^Rf^L) 
+ 9lRieL'^eR)0^Lf^R) + 9lLieL'^eR){'i^Rf^L) 

+ 9lR{eRY'^eR){'i^Rlt,f^R) + 9lLieLY'^eL){'i^Llt,fiL) (109) 
+ 9IR{eLY^eL)0^Rl^^^^R) + 9ILi'^LY'^eL){'i^Rl^,f^R) 

9^ - — 9'^ - — 1 

+ -Y^ieRa'""UeL)i'l^RCri,^fiL) + -^{eL(yi,uVeR){^L(^i,vf^R) + H.c. 

The Left-Right Symmetric Model (the LR Model) is an extension of the SM to restore the 
parity symmetry in the original Lagrangian. The LR Model is based on the gauge group 
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of SU{3) ® SU{2)l ® SU{2)r ® U{1)b-l- The charge of U{1) in the LR model has a clear 
meaning as the difference between the baryon number B and the lepton number L in contrast 
with the mysterious hypercharge Y in the SM. Similarly to SU(2)l doublet in the SM, the 
right-handed fermions compose doublet of SU{2)ii in the LR Model. Therefore, the right- 
handed neutrinos Uji are introduced naturally as SU(2)ji partners of right-handed charged 
leptons. After the spontaneous breaking of SU{2)r(^U{1)b~l to U{1)y, the hypercharge is 
given by 

Y/2 = hR + {B-L)/2. (110) 
Since electric charge is connected by 

Q = l3L + Y/2, (111) 

f lllOP implies the charge quantization, which is one of great achievements of {sSj. 

Since we require the parity symmetry to the theory, the gauge coupling of SU{2)ji must 
be the same as the one of SU{2)l: §2 = g2L = 92R- The Higgs field which gives the Yukawa 
terms is a complex bidoublet of SU(2)l ® SU(2)r with B — L = 0. The bidoublet field can 
be expressed as 

/ (A? 

(112) 

which transform as = Ul^UJ^ under SU{2)l and SU{2)r. 

(k o\ 

($) = 1 (113) 





with K ^ hi' gives rise to the breaking of L-R symmetry. However, $ is neutral (B-L=0) and 
U{1)b-l is not broken. So we need another Higgs. Usually, two complex triplet fields (A/, 
for SU(2)l and A/j for SU(2)r) with B — L = 2 are also introduced to generate Majorana 
neutrino masses (see also the Higgs Triplet Model in Sect. lIVEp . 

Ai = (3,l,2) A^ = (l,3,2) (114) 

Then, the gauge symmetry breaking proceeds as follows: first A^j acquires vev vr, leading 
to SU{2)l X U{1)y with flllOp . which furthermore breaks to U{1)q by the vev of 
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FIG. 13: One of diagrams of gauge boson contributions to the electron EDM in LR symmetric 
model. Gauge bosons contribute to quark EDM also. 



The triplet Yukawa coupling for Ar must be equal to the coupling for Al because of the 
parity symmetry which is spontaneously broken by v^. Thus we have 

Vr > k, k' > vl (115) 

Thus the two Higgs doublets model (not of all but its measure part) and Higgs triplet 
model in the previous subsections are combined together in left-right model. 
Their vevs vr and are different to each other and from k, and k'. 
Figure [13] shows one of diagrams which contribute to the EDM of the electron. In a 



simple case where there is only one flavor of leptons, the electron EDM is estimated 57| as 

2 



_27 \lm{mD)\ c_ [ niR 



8.2 X 10^^' ' ' e cm for > 1, 

141 <<; \^wj (116) 



-26 |Im(mB)| ^ _ ( niR 



3.3 X 10"^'' ' ' " e cm for ^ < 1, 
MeV \mw J 

where mn denotes the Dirac mass of the neutrino {mDiTZvR) and m^, mw = 80 GeV are 
the masses of heavy right-handed neutrino and the light W boson, respectively. 
The contributions of Higgs bosons in the LR model are not significant 



Another important contribution of Wr in CP violation may be the neutrinoless double 



beta decay [59 1 
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G. Fourth Family Model 



We set the quarks of the fourth family [60|] as 



it',br. (117) 

The mixing angles and CP violating phases for N families are given by 

9 . ^ N(N-l) (N-l)(N-2) , , 

iV^ - (2iV - 1) = '-^^ '-^ (118) 

where the first term is mixing angles and the second CP phases (see Appendix D for Majo- 
rana fermion case). 



If we consider 4-generation SM (SM4) 26[, we can construct new Jarlskog parameter in 
place of f l53|) 

^(234) = - m^^){ml - rr?^{m\ - ml){ml - ml){ml - ml){ml - ml)J(^2M). (119) 

If we take heavy quark masses t' and h' in the range of 300 to 600 GeV, ^(234)/7'^vf 
of order UB/n^- 

\Vud? + \Vus? + \Vub? = 0.9999(10). (120) 

Thus SM4 enhances CP violation and therefore the EDM also. 
We may consider inside the loop only t, t', b' heavy and identify as 

c = u = U, d = s = b = D. (121) 

Then two loop diagram of Fig. 2 vanishes for u quark but survives for d,s quarks, giving 

d, ^ X'^^T^.Gjm,^ ^ 3 X lO-^^ecm, (122) 
Air Air IQir^ 

where A is the Wolfenstein parameter, A = \Vus\ = 0.22. This is only two orders of magnitude 
larger than the SM result of section fllll A3p . However, if we consider the chromoelectric 



dipole moment of the s quark. 



~ .as aw 1 ^ 

- X'^T^.G,ms% (123) 
All An IGvr^ mi 



Using the estimation of the relation between and ds by 6l|, we have 

dN ^ ~ds ^ 5 X 10~^° e cm. (124) 
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H. Extra Dimensions 



The motivations for extra dimensions are diverse for both SUSY and non-SUSY. There 
are many SUSY breaking scenarios. The extra dimension makes the geometrical SUSY 
breaking possible. The gauge supermultiplets propagate in the bulk, we get gaugino mass 

" wk- ''''' 

which is called gaugino mediation. 

Even if the theory itself is CP invariant, it may be violated by extending the theory to 
extra dimensions. This is because the compactification of the extra dimensions does not 
respect the symmetry in general. The CP phases come either from the boundary condition 



of extra dimensions (One of Scherk-Scwarz mechanisms 6^) or vev of fifth gauge field 



(Hosotani mechanism 63|) 6J| 65 1. 



L = iilj-i^{dN -ieA^)^ - M^^. (126) 
Inclusion of a torsion term results in a nonminimal term 

K^aFip (127) 
or the Fermion mass term via Hosotani mechanism 

^(M + ^75^4)^^ (128) 

with 

Xi= I dyA^. (129) 



Here y is the coordinate of the extra dimension. The rotation of mass term gives rise to 

K'^aF'y^ip. (130) 



The concrete constraints from the observation are given, for instance 



66| 



4 14 
dn = -dd - -du = -dd, (131) 

simply because they did not consider up quark and 

d{KK) ~ -2.3 X IQ-'^^iRmwY [e cm]. (132) 

Since ^d{KK) must be less than the experimental upper limit, we have 

4 > S3mw ^ 2.6 [TeV]. (133) 
R 
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V. THE EDMS OF ATOMS 



The origin of the difficulties of the measurement of electron EDM is due to the absence 
of resonance unlike the neutron. A possible way is to perform the resonance experiment on 
neutral atom and to interpret the result in terms of electron EDM or hadron EDM. These 
object has very tiny values and let's consider the effect linear in the EDM. In the subsequent 
atomic and molecular experiments we treat an internal electic field Ej^^ induced by atom 
or molecule as well as an external elctric field E. This E induces the EDM er^ with the 
intrinsic = Xli^e- Thus the total Hamiltonian is a sum of unperturbative P,T-even 

term, 

Ho = ^c(x-p^ + /3,mc^ + Vnuciiri) + J] Vc{r,j), (134) 

i i<j 

and T,P-odd term 

HPTV = -E^e ■ EL - 5^d^ E - ej^r, ■ E. (135) 

i i i 

The last term of Hq is a two-body interaction and can not be solved exactly. 
The first and third terms are P-odd and the second P-even. So the first and second order 
energy shift are given by 

El = -J2M^e\mo)-^ (136) 



and 



e;2 _ f ("^o|de • Eint\no){no\er' ■ E|mo) (mo|er^ ■ E|no)(no|d^ ■ Eint\mo) ) 



n^m i 



(137) 

Here |mo) is an eigen state of Hq, it should be remarked that, as will be shown in f l330p . 
EDM appears as the coefficient of the energy shift linear on the external electric field. So 

E^ = El + El = -d' ■ E, (138) 

where 

(mold* ■ Eint\no) {no\er'\mo) 



EO - E^ 

m n 

i 
e 

El - m 



(mo|er*|no)(no|d*e ■ Ei„t|mo) , 
7^ T^n r • \^^^) 



m n 
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However this d' vanishes as follows. 



(mo|er* ■ E|no)(no|d^ ■ Ei„t|mo) = ie{mo\r' ■ E|no)(no|d* ■ [p*, Ho]\mo) 

= ze(mo|r^ ■ E\no){no\dl ■ p|mo)(^°, - (140) 

Using the communtation relation = iMjj the second term of fll39p cancells with the 



first term. This is the famous Sciff theorem 



67|. Since the expectation value of 



S ■ E,„i = [S ■ V, Ho] (141) 
does not contribute to a linear Stark effect, the residual EDM is becomes 

Vedm = -4(/3-l)5]-Ei„i (142) 

and the residual energy shift is 

AE = -de{mo\{(3 - 1)S ■ E|mo) - 24 2^ _ 

= -d{atom) ■ E (143) 



The first term has no enhancement factor unlike the the second term and is much smaller 
than the second term, and 

diatom) = -24 E E ("^ol^>°)<^°l(/^-^)^^-E-^l"^o) (144) 

datom has the large value when these states are almost degenerate. However, this enhance- 
ment is reflected in quite different ways in paramagnetic atoms and diamagnetic atoms. 
Though 01431) itself is rather universal, Hp^v is variant. One example is P,T-odd Nucleon- 
electron interaction like (see Appendix C for the detail) 

+ iGs'NNL-f^L + iGp'N-i^NLL. (145) 

There are the other CP violating effective interactions (see the last part of section IIB). 
Another important interaction is due to Schiff moment. The origin of Schiff moment itself 
is not unique. 

As we mentioned, in the nonrelativistic Hamiltonian for a system of particles of finite 
size, there is no interaction energy of first order in the EDM if there is no misalignment of 



charge and moment distribution. Schiff also indicated in 
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67l | how this theorem is violated 



by relativistic (Breit equation O ((f/c)^)) and the misalignment (the Schiff moment), where 
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indicated that 



V is the velocity of electron or nucleon. Prior to this discovery, Salpeter 
radiative corrections of O {{v/cY) enhances Hydrogen EDM. 

Sandars pointed out that relativistic effect of electron EDM in heavy alkali atom gives 
large atomic EDM 69|. 

Let us proceed to discuss in more detail for Hydrogen-like atom. The EDM has, by 
definition, odd parity and naively vanishes between the states with same parity. 

For nonrelativistic case, its energy levels are 



E 



2n2 



(146) 



Here n is the principal quantum number, and this energy is degenerate n^-ply, = 0" ^(2/+ 
1) = ^2 (see Eq.f lT59|) l 

If we consider the relativistic effects (spin effects) the degenerate energy levels are split 



into n fine-structure components at different j [70|. Let us obtain the relativistic terms 
w.r.t. 0{v/c) (see Appendix E for relativistic expansion and Appendix F for nonrelativistic 
approximation for more detail). 



At the first order of v / c, we obtain the Pauli equation 



P 



-a B 



2m \ c J 2mc 
In further approximation of O {{v/c) — 2), we assume B = {i.e. A = 0), we get 



2m 



6$ 



eh 



If E is centrally symmetric, 



E = - 



Am^c^ 
r (i$ 



cr ■ [E X p] 



eh^ 



■im^c^ 



V-E. 



r dr 



The spin-orbit interaction (the fourth term) becomes 



eh 



si 



Am^c^r 



cr ■ [r X p 



rf$ _ h^ dU 
dr 2m'^c'^r dr 



Is. 



For many electron case of atomic number Z, 

Vsl = ^ aJaSa, 

where 



an 



h^ dUjrq) 
2m'^c^ra dra 



(147) 



:i48) 



(149) 



(150) 



(151) 



(152) 
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and, therefore, 



ttB 

2 \ 2 4 

^ me 



a 



he 



For given total L and S, the averaged Vsl is 



Vsl = AS-L, 
L ■ S = hjiJ + 1) - L(L + 1) - S{S + 1)]. 



(153) 
(154) 

(155) 
(156) 



Since the value of L and S are same for a multiplet, energy splitting is given by the Lande's 
interval rule. 



AJ. 



(157) 



Then we obtain 



(2Si/2, 2^1/2) , 2p3/2 

(3si/2, 3^1/2) , (3^3/2,34/2), 34/2- 



(158) 



The remaining degeneracy is removed by the hyper fine-structure components caused by the 
radiative correction (Lamb shift Zlj). So using this hyper fine splitting, we can obtain large 
atomic EDM 



To estimate atomic EDM we need two informations. One is that of atomic wave functions 
and another is that of P (or T) violating interactions. 

Atomic EDM are due to those of constituents, electrons and nucleons. For electron 
EDM it is very important that the atom has an unpaired electron, and electron EDM is 



proportional to 



69j (For the review, see [10|[1]J].). If there is no unpaired electron 
(diamagnetic atom), we can measure quark (or hadron) EDM. For proton EDM, nucleus 
has an unpaired proton. In this case polarized molecule takes an important role 72|. 

We are dealing with many electrons system and the electron wave functions are not in 
general exact. In this case the expectation values of the EDM depends on the representation. 



< 6|r|a >= Tfta 



1 



Eh — En 



< b\HQr - rHo\a >, 



(159) 
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where Hq = — + V. Inserting this into fll59p . we obtain 

na = -ET^—ET < V'r - rVV >= —pba (160) 

= -^{VV)ba. (161) 

These three representations are of course equivalent. However, if we use the approximate 
wave functions, these values are different in General. So we must be careful what is the 
oHgins of d,sa.epanc,es. due to d.ffe.e.t app.ox„„a«ons o. to .ep.ese.tat.ous fl- 

A. Relativistic Effects 

The relativistic equation of atom with CP violating interaction (^ term) is 



I *e \ .mc 



« = ^^^l5l^l.F^.u. (162) 



Electron EDM breaks the CP invariance and CP violating energy equation is 

{Ho + H')u = Eu. (163) 

Here Hq is the Hamiltonian of the original single electron Dirac equation in the external 
field, 

Hq = m(3c^ + a ■ (cp - eA) + e0 (164) 
which leads to f ll34p in the static limit and H' is CP violating interaction Hamiltonian of 

6 



the right-hand side of 



68 



where S is defined by fl27|) and 

/ rr\ 

(166) 

^ is dimensionless constant which measures the EDM in units of the Bohr magneton. Thus 
^ -C 1 implies that the EDM is small compared with e times the Compton wavelength. The 




This term only contributes to the intrinsic EDM and we denote H' , distinguishing it from HpTv of (|135p . 
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last approximation in fll65p comes from the relativistic suppression due to the mixing of 
the upper half with the lower one. It should be remarked that cr ■ E is T-odd (T is time 
reversal operator) but a ■ B (the suppressed term) T-even. 



We consider a hydrogen-like atom with charge Z, where E = ^e. 



H' = —^Zar "^Sr with Sr = cr ■ er/2 



(167) 



in atomic units e = m = h = 1. The Hamiltonian of the above single e 
mean field approximation is generalized for many electrons system as 

2 



ectron system in the 



69| 



B 



•jk 



(168) 



Here suffix indicates the quantity of i'th electron and Bjk is the relativistic corrections. 
Taking the retard potential into consideration, Lagrangian can be described up to order 
0((^)^) as 



h Bjk = — 



rjk 



'jk 



(169) 



further higher order > O ((f )^) corrections come from photon emission (Breit interaction) 
and 



mc 



?fcPjr)(rifcPfc) ^ 



-{a-j + 2(Tk)[rjkPj] + {(Tk + 2(Tj)[rjkPk]) 



(170) 



1 / eh\'^ ( (TjCTk ((Tjrjk){crkrjk) Stt 

— ) I ^ - ^^^5 — - —(T,(TkS{T,k) 



'■jk 'jk 

The second line corresponds to spin-orbit interaction and the third spin-spin interaction. 
If we incorporate the spin of nucleus, the degeneracy of J is split (hyperfine structure). 



Vij = ai ■ J, 



(171) 



where i and J are the spin of nucleus and total angular momentum of electron envelope, re- 
spectively. However, in this hyperfine splitting dominant contribution comes from magnetic 
dipole and electric quadrupole and does not play important role in the EDM. 
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Thus the hnear Stark appears as relativistic effects in duphcate meanings, i.e. 1-/3 and 
Bjk components. Many particle interaction effects are due to this relativistic effect as well 



as due to the other nonrelativistic excitation effects. 

We proceed to the detailed calculation of single electron case f ll67p 
commutes with and M^, but has odd parity, and 

<l = J± ^,J,m\sr\l = J T ^,J,m >= -. 
Let the radial part of m = rxni- Then it satisfies 

{ [2fi(P - S/dr^] + /(/ + l)r-2 - 2^Eni] Xm = 0, 



The operator Sr 



(172) 



(173) 



where fi is the reduced mass of electron and set equal to unity in the subsequent equations. 
It follows from ffT72|) that 



<n,l+,j,m\H'\n',L,j,m>= -^Za{2j + 1) \En - E'J J drXni+Xnri. 



;i74) 



with /-t = j ± i. Therefore, naively the first-order perturbation vanishes. The exception is 
discussed in section 5.2. The second-order perturbation energy is obtained by use of (11740 



AK,z±,, = {^Za/{2j + 1))2 J drxn,i^ 



1 rf2 Z l/^(/^ + l) 
J^n + - -rr; + 



Using 



we obtain 



2 dr^ r 2 

^2 



Xnl± ■ 



n3(/ + l/2)' 



Ai?„,;^,, = ±Z'en-%2j + l)-i(/± + ^)-WRy 



with Ry=rn^ = 13.6eV. 



B. Peculiar Property of Paramagnetic Atom 



The atomic enhancement factor defined by 



(175) 

(176) 
(177) 



is given by \L 



(J + l)a|7(47' - l){NoNraf/\E„, - E,) 



:i78) 



(179) 
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for alkali atom. Here the sum is taken over the excited state m, and Nq, N„i are effective 
principal quantum number defined in fllQOp . 7 = ^ ( j + 1/2)^ — Z^a^ and rom is electric 
dipole radial integral, 

POO 

rni,n'i' =< nj\r\n',l - 1 >= Vl Rn\i-iRn,ir^ dr (180) 

^0 

in units of = ^|^:Bohr radius. We will derive f ll79p (see Eq. fll93p ). We start with the 
general relativistic arguments. For Diamagnetic atoms the dominant contribution to atomic 
EDM comes from that of nucleus, which will be discussed later. 

In hydrogen-like atom, states of different angular momentum / with fixed principal number 
n are degenerate in nonrelativistic approximation. The eigen functions with external field 
are the superposition of the field-free function with different 1-values, which gives the linear 
Stark effect. Let us write the Dirac spinor in the form 



u± = r ^ {x2±{r)Vji±, - iXi±{r)VjiTy 



Xi satisfy the following equations. 



dxi 
dr 

dX2 
dr 



_Xi 
r 

X2 





mc 




E \ 


z' 

— a— 










mc 




E \ 


z' 









X2 



where 



= T(j + ^) for j = / ± ^. 



Using (11651) and (I172p . we obtain 

{nJJ+,m\H'\n,j,l^,m) = --^Za J rfrr"^(x2+X2- - Xi+Xi-)- 
X2± and xi± are related to each other as (I182p and we obtain 

POO 

4 {h\H'\l_) = iZa' / dr{D^x+)r~\D^X-), 
Jo 



where 



d + 1/2 
dr r 







The exact Dirac wave functions with given n,l,j are 

X2 ^ Vr(27 + n, + l) / 1 + 6 /2Z_^ '"^ 
r V{2-i + l)^. ^ 4.N{N - K)\NaB 

2Zr 



2Zr 



Na 



B 



7-1 



X 



-UrF 



1,27 + 1, 
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{N -k)F -r2„27 + l. 



2Zr 



Na 



B 



181 



;i82) 



;i83) 



;i84) 



(185) 



(186) 



(187) 



and 



Xi 
r 



1 -e 



2Z 



v/r(27 + + 1) 

r(27 + l)v^ \^N{N-k) \NaB 

2Zr 



3/2 



Zr 



2Zr 



Na 



B 



7-1 



X 



rirF -rir + 1,27 + 1, 



Na 



B 



{N - k)F [-nr,2-f + l, 



2Zr 



Na 



B 



(188) 



Here F is the confluent hypergeometric function and radial quantum number, tlie number 
of nodes of radial part of the wave function, 

aZe 



7, n = rir + \k\ 



with e = and 



Here we have used 



N 



n 



^ -2nr{\K\-^K^ -a^Z^). 



2 '72 

a Z 



;i89) 



(190) 



(191) 



and normahzed Ar = as in the hypergeometric functions. N is called apparent principal 



quantum number and 



mZ'^a^ 
2iV2 



E. 



nl 
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Using these forms, fll84p finally reads, for instance 

^Z^ai-f - 1 



(2Si/2|if'|2pi/2> 



27(27 -1)(7 + 1)(27+1) 



(192) 



(193) 



where 7 takes the value y/1 — Z'^a^ in this case. For small Za, it is reduced to 



(2Si/2|i/'|2j9i/2> = 



For more general case 



<j,/+|i/'|j,L >= 



8\/3 



4(Za)= 



(194) 



7(472 -l)(iV+Ar_)3/2 
For heavy alkali atom, for instance, cesium, {9] 

16 Z3a2r(6s,6pi/2) Ry 



Ry. 



(195) 



K{Cs) = d{Cs)/de 



3 a|7(47' - l){NsNpy^/^ ^(6^1/2) - E{6s) 
The radial integral is experimentally known 77 1 



118. 



(196) 



/■oo 

r(6s,6pi/2) = / drr^RQo{r)RQi{r) = S.Sa^j. 
Jo 



(197) 



41 



fll96p should be checked with the experimental result [24 1. 

For Francium (Z=87, 7s— I'Pxj'i)^ K{Fr) is estimated as 873. However, as was stated in 



10|, fll79p is not applicable for atoms with complex configurations and requires electrons 



78| and K{Fr) is modified 



correlations. Such calculations are performed in, for instance, 
to 895. 

There are some discrepancies on the estimate of enhancement factor of K of Thallium 
[Xe]4/i45(ii°6s26pi 



sumptions. 



79| 



8l| . The discrepancy seems to come from the starting as- 



considered that Thallium has three va. 



ence electrons, 6s^6p^, whereas 81 1 



considered it having one valence electron. If we adopt |80| 

d(205T/) = -582(20)4 



or if we take 81 1 



-4664 ^ 4 < 1.6 X 10~^^ 



e cm. 



(198) 



(199) 



82 1 which asserts 



In preparing this revised version, an interesting paper has just appeared 
that this discrepancy disappears, converging to K = —573. 

Xe has closed electron shell of 5s^ 5p® but we may one elctron of 5p state excited to 
5p^ 6s^, which resembles with that of Cs, [Xe]6s^, whose enhancement factor was estimated 
to K{^^^Cs) = 114 3 or 120.53 84]. 

As for ^^^Xe, the lowest excited state with a 6s electron has a enhancement value 



K{^^^Xe*) = 120 [75||85| or 111 [86|, and 



4 < 3.2 X 10 



-24, 



e cm. 



(200) 



Using these results, Ellis et al. considered the maximal EDMs of nuclei 



C. Chiral Condensate 



Before discussing diamagnetic atom, we will briefiy resume QCD chiral dynamics [88|. 
This is because hadronic matrix elements are described in terms of quark condensates by 
using operator product expansion. 

Let us begin with the following effective action (see Appendix G for the implication of 
the effective action). 

L = q{YD, + M)q - ^G,,G^\ (201) 
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where 

D,, = d^ + z^7.A;:A^ (202) 

This action is invariant under SU{3)l x SU{3)r transformations in the hmit of = = 
rris = 0. That is, 

Q^g = e-"^'^g, Qj,q = e'^''^''^'q, (203) 
where u, d, s quarks constitute SU(3) group, 

q = {u,d,sf (204) 

and we have the following conserved currents 

f^R = QL,R>^''rQL,R- (205) 

Here A" are the Gell-Mann's 3x3 matrices and ql (qr) are left-handed (right-handed) part 
of q. 

= jT+Jr', (206) 
= jT-Jr'- (207) 

So we have the conserved currents and conserved charges Qa and Qsa- They satisfy the 
algebras 

[Qa, Qb] = ifabcQc, [Qba, Qb] = ifabcQbc, [Qba, Qbb] = ifabcQc- (208) 

However, this group is not exact and they are spontaneously broken to 

Qa|0)=0, g5a|0)^ 0. (209) 

Thus there appear 8 pseudo Nambu-Goldstone bosons. Pseudo implies that the original 
chiral symmetry (Qs transformation) is not exact. It is broken by 

HsB = rriuUU + m^dd + rrisSS. (210) 

This can be rewritten as 

HsB = {rriu + rrid + ms)(uu + dd + ss)/3 

+ {rriu- md){uu -dd)/2 (211) 
+ {2ms — iTT'u — iTT'd) (2ss — uu — dd) 1 6. 
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FIG. 14: Chiral anomaly induces vr" — )■ 77 via ifonFF interaction. 



Here the first line is an SU(3) invariant, the second breaks isospin SU(2), the third responsible 
for the SU(3) symmetry. 

M2 = (m„ + md)B + 0{m\ In m,), (212) 
M|.± = (m„ + m,)5 + 0(mj In m,), (213) 
M|.o = (md + m,)5 + 0(mj In m^). (214) 

Here B = — |(0|gg|0) with pion decay constant /^^ = 133MeV, and we have used the chiral 
limit 

U = Ik. (0|mm|0) = (0|drf|0) = (0|ss|0). (215) 
Adler-Bell-Jackiw axial vector singlet current anomaly 



89| and its non-Abelian version is 



d,]'^ = 2z rn,q^^q + ^ (^FF + G^G'') . (216) 



N 

q=u,d. 

with the number of flavour Np. Since 



87r2 
with 



QaQa ^ 2Nfd>'K^ (217) 



1 / 1 

cabc Apb Acre 



K, = ^e^^p^A-- (^aM™ + -gr^AP'A'^'J . (218) 
In the limit of rrtu = 771^ = rris = 0, the axial current j^^ gets conserved again by replacing 



j^^ as 



Jsm = - 2NfK^. (219) 



Thus quark condensate is essential for chiral symmetry breaking. Nevertheless, anomalous 
term is crucial to the presence of 7r° — j- 27. The neutral axial vector current gives the 
modified PCAC relation. 



5% = Uraln' + —FF. (220) 
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We will come back to this problem in the next subsection. U(l) problem concerning QCD 
condensate is discussed in Appendix H. 



D. Peculiar Property of Diamagnetic Atom 



Diamagnetic atom has no unpaired electron, and the main contribution of atomic EDM 



comes from the misalignment between charge and the EDM distribution of nuc 
hadronic part of the atomic EDM manifests itself through the Schiff moment 



eus. Thus 



67|(90| 



Hatom -^eZectron ~l~ H nucleus ~l~ ^ ^ (^'^(^0 ' ^nucleus ' E, 



(221) 



where are the i'th electron coordinates and dnucieus is the nuclear EDM. E means the 
external E. It should be remarked for diamagnetic atoms that the Stark term due to 
electron' intrinsic spin term (the first term of (11351) ) is replaced by that of nucleon. $(r) is 
the nuclear electrostatic potential given by 



<l>(r) 

where p(x) is the charge density of nucleus. 
Here it is important to notice 91 1 

" z 



(222) 



— t 



'-atom 



4 = 1 



:^V,<l>(l 



ZeE, 



(223) 



i=l 



where pj are the momentum of atomic electrons, and the first term is the average electric 
field induced by atomic electrons. The expectation value of this commutator in the energy 
eigenstate vanishes and we may add 

1 ^ 

V = dnucieus ■ E - ^ ^ dnucieus " ^i^i^i) (224) 



1=1 



to Hatom. as far as we consider the expectation value. This implies we may change 



-d 



nucleus 



E ^ -(d 



nucleus 



-id 



nucleus 



))-E. 



(225) 



So the expectation value is zero. This is another statement of the Schiff theorem. From the 
first term of fl223p . we should consider the interaction of atomic electrons with the nucleus, 

1 



$fr) - ^ {dnucieus) " V$(r) 



(226) 
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as the screened electrostatic potential. Therefore, the atomic EDM reads 
0|eEf r^lri) (rije^f ($(ri) - (d^e^s) ■ V$(ri)) |0 



'■atom 

n 



^ + (227) 



Using the charge distributions 

j p{x)d^X = Z\e\, J Xp{x)d^X = {dnudeus) 



y a;V(x)rf^x = Z|e| (x^)^^ , J {xkXk' - ^6kk'x'^)pix)£x = Z\e\ {Qkk') etc (228) 



N 



e<l>(r) - ^ (dnudeus) ■ V$(r) 



Ojv ) = -— +47reS- V5(r) + ... (229) 



Here ... indicates electric octupole and higher pole contributions, and S is the famous Schiff 



moment 



92| . (The detailed derivation is given in Appendix B 

z 

10 



(230) 



The (Qkk') vanishes for ^^^Hg, ^^sxe, 225Ra. 

There is another Schiff moment S""'^' due to the misalignment between the charge distri- 
bution and the EDM distribution of nucleus, whose derivation is given in Appendix I. 

Corresponding to these situations, we should consider fl230p more generally 

S = 1^ E E ((^^ + P^^" - \ (^')^'^) + P^- (231) 

Here r^r is a N'th nucleon position and pj is the position of the ith charge inside the N'th 
nucleon, and 

^ Cj = Cat, ^ Cipj = djv. (232) 

i i 

Retaining the term up to linear term of p, we have 

S = S^^ + S"""^ (233) 

where S'^'^ is given in f l230p and 



N N ^ ' 

Usually S""'^' is considered to be small compared with S'^^ . The mean value of S'^'* is 
nonzero only in the presence of P- and T=odd nucleon-nucleon interactions. 
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(234) 



FIG. 15: One 5^*^ coupling induces effective CP-odd NN interaction, which give rise to S"^. 



N 



FIG. 16: One g^''"'^ coupling induces d^v, which gives rise to S 



nucl 



For the arguments of hadronic EDM, we must represent hadronic CP violating interac- 
tions from those of fl39|) and fHOjl . This will be discussed in the last part of this subsection 
(see (12801) ). They are described as 

C.NN = g^j^Nr'^Nn'' + Q^^LnNtt' + (^^^(iVr^iVTr'' - SNr^Nn'). (235) 

Here g^j^j^f {i = 0,1,2) are CP odd coupling constants, whereas we denote the CP even 
strong ttNN coupling constant as Gt^nn{= 13.5). The Schiff moment due to this coupling 
is calculated as follows. f l235p gives rise to both S"'' and P- and T-odd NN potential 

has the form via FigHJ] and its effective potential is given by 

W{r^ - r,) = [^(0)(^^ . r,) - ^(r^ + r^) + g^'\3r:n^ - r. ■ r,)] {a^ - a,) 

Tri'^a -n)i(Ta + a-b)} ■ {ra-n) . ^ 1 + . r . 236 

Here we have suppressed the subscript ttNN in g^^\ dj is via Fig {16] and 

In 

Given T and P-odd perturbation, let us calculate Schiff moment 



p™^^-^ (o)_ (2)) (237) 



m 



H = Hq + Hres- (238) 
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Here 

i^o = T + K)o + Vii (239) 
is unperturbative one-particle Hamiltonian and exactly solvable and 

Hres = W + V22 + + + V^04 + ^40- (240) 



W is the pseudoscalar interaction fl236p and V the Skyrme interaction [9J]. Subscripts {ij) 
refer to the final and initial numbers of quasiparticles. 

Let us assume that in the O'th order approximation, the state is = |«), and define Q 

by 

Q^Y.\^)m. (241) 

Then perturbed wave function is given by 

*a = f 1 + ^ ^ Hres + ^ ^ Hres ^^ Hres + $ a- (242) 

This is the Brillouin-Wigner expansion and is the single quasiparticle energy of the valence 
nucleon. 

So in the first order perturbation of S*^, we obtain 

{'^a\S'\'^a) = iV-'($a| [1 + Hres { + -] [l + Hres + •.•] |$a). (243) 



The first-order quasiparticle (Goldstone) diagram is given in FigJTT] [95 1. 
Here the Goldstone diagram implies that 




(244) 



liger order quasiparticle calculations need some elaborate code and should be referred to 



97[, and we simply list the final results 

i^alS'l'^a) = S=iao + b)G^NMg^^^ + a.G^N^g^''^ + (aa - b)G^NN9^'\ (245) 

where the coefficients a, specify S'^^ and b does 5'""'^' defined by f l233p . The numerical results 
of ai and b for ^^^Hg are given in Table [11 

In this Table, the first two papers considered that EDM of nucleons d^v is independent 
on C-„NN, whereas Ban considered djv is related as (123 7p . The second and third papers 
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(i) 



(iv) 




(ii) 



-X 5i\ 



/wv^X Wii W20 X'wW 



(v) 
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(iii) 



W31 



(vi) 



X >S2{ 



X Sn' 



95( 1 . The broken 



FIG. 17: First-order quasiparticle diagrams contributing to the Schiff moment 
hne represents the action of the Schiff operator, the zig-zag Une represents the P- and T-odd 
interaction. The looped line in higher order diagram represents a generic Skyrme interaction. 



TABLE I: Calculated coefficients Cj and b for ^^^Hg. The units are e fm"'. The last two references 
include the Skyrme interaction SkO'. Five results of Ban et.al. are due to Hartee-Fock and 
Hartree-Fock-Bogoliubov approximations. SLy4, SIII et al. indicate several Skyrme interactions. 



Dmitriev-Sen'kov 2003 [96] -0.0004 -0.055 0.009 - 

de Jesus-Engels (averaged) [95] 0.007 0.071 0.018 - 
Ban et al [97] 

SLy4(HF) -0.013 0.006 0.022 -0.003 

SIII(HF) -0.012 -0.005 0.016 -0.004 

SV(HF) -0.009 0.0001 0.016 -0.002 

SLy(HFB) -0.013 0.006 0.024 -0.007 

SkM*(HFB) -0.041 0.027 0.069 -0.013 

incorporated a collective modes based on the different approximation methods for nuclear 
structures but the results are still divergent. Thus it is difficult task to precisely estimate 
the EDM of diamagnetic atom and to induce nucleon or quark EDMS due to lack of precise 
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theory of nuclear structure. Taking this point in mind, let us consider some cases. For ^^^Hg, 
Numerical calculation is 



rf(i99Hg) = -2.8 X 10 



-17 



S 



e cm. 



(246) 



The value of of the Schiff moment of d{^^^Hg) can be presented as a sum of proton and 
neutron EDMs |99| 

(247) 



S Spdp -\- Sjidfi 



with Sp = 0.20 ± 0.02 fm^ and s„ = 1.895 ± 0.035 fm 



Combining the experimental value 



100| (see more up-to-date data in 23|) 



dC^mg) < 2.1 X lO^^^e cm 



(248) 



with f l247p . we obtain 



\dp\ < 3.8 X 10"^^ e cm < 4.0 x 10~^^ e cm. 



For ^^^Xe case, numerical calculation is 101 1 



dC^^Xe) = 0.38 X 10 



-17 



s 



e fm^ 



e cm 



The measurement is 102 1 



(249) 



(250) 



dC^^Xe) = (-0.3 ± 1.1) X 10-^*^6 cm. 



26, 



From (12511) value, |l03j obtained 



(251) 



\dp\ < 4 X 10 ^^e cm \dn\ < 1 x 10 ^^e cm. 



(252) 



General arguments of CP violating four Fermi coupling are given at Appendices C and D. 



1. cEDM and parity odd nuclear interaction 

In this subsection we give a very short review of chiral symmetry and its breaking in 
strong interactions since it has many problems. Let us start with the conserved axial- vector 
current (CAC) hypothesis jl04 |. 

d,jU^) = 0. (253) 
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FIG. 18: Nambu's interpretation of pion dominance. 



Of course GAG requires = = and is not realistic. However, it makes clear to 
understand how to break chiral symmetry. f l253p leads us to 



{N'\J^,\N) 



m 



N 



EE' 



Nambu aserted 



105| that 1/g^ in the second term of (12541) should be interpreted as 



1 

q2 



lim 



1 



This corresponds to the diagram (Fig. (!T8|) ), which can be written as 

Comparing (I256p with (r254p where (I255p is inserted, we get 

9v 



This is the Goldberger-Treiman's relation jl06 |. Here use has been made of 



2uj 



(254) 



(255) 



(256) 



(257) 



(258) 



{N'\d,jr\N) 



—I 



m 



N 



EE 



-Mt) 



-2mNU 75 — u — 2m 



N- 



q'- — 



Substituting the equation of motion of tt, 

(□ + m^)7r-=r, 



into f l259p . we obtain 



(A^'|7r"|iV) 



{N'\3-\N) 



mt 



EE' 



gZ _ 



(259) 



(260) 



(261) 
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Assuming the matrix elements vary little between t = and t = m^, we obtain 

{N'\d,jrmN) ^ ^-^ml{N'\n-\N) 

QV ^ttNN 

= Uml{N'\7i''\N). (262) 

Thus we obtain the PCAC condition 

d,jr = Umln'^ (263) 

and (12201) with chiral anomaly. 

Next, we proceed to discuss the path from the presence of the strong EDM of dimension 
5 (cEDM and 9 term) to the effective CP-odd g^^. 

Let us write hadronic CP-violating operators like fl39l) and fHOj) etc. as O and consider 
the following two points correlation function ''. 

q,M, ^q,J d'xe-''^^N\T 0^^)0(0)) \N'). (264) 

From the definition of time ordered product, the right-handed side is rewritten 



q,M, = -tj d^'xe-'^^HNlT {d,j^^{x)O{0))\N') 

-z6{xo){N\[j^„ Om\N')}. (265) 



Using (12631) and LSZ reduction formula lOSi ] in g — >■ limit, we obtain 



q,M, = -Un^N\OmN') - z(iV|[Q^(0), O(0)]|iV') (266) 

or equivalently 

limy2^(vr"iV|O(0)|iV') = -f (iV| [Q^(0), 0(0)] |A^') 

In 

-lim| J dV^N\T{J^^{x)O{0))\N'). (267) 
Substituting the concrete form of 0(0) as (HOj) into the above equation and using 

[Q^(0),g(0)] = zt"75g(0) (268) 



^ The general arguments on the operator expansion of T product of two currents are given in [10/ 
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with the generators of group of flavour t", we obtain 109l |. 110| . 111 | 



RHS of dMD = —{N\du{uGau - mluu) - dd{dGad - mldd)\N) 

Jit 



+ — 



with 



niu nid vfis 



niumdms 



{N\uu - dd\N) 
rriumd 2 (OIMCo-^IO) 



, mo 



mQ is estimated as 



from QCD sum rule 



ml ^ O.SGeV^ 



(269) 



(270) 



(271) 



Peccei-Quinn mechanism 



112 L; Here quantum corrections are also included. If we us e th e 



113| . the second term of fl269p vanishes in the following way |114 |. 



L = —aGG, 
Stt 



(272) 



where a is axion field and GG = G^^^^G^^p. When there exists cEDM, axion potential 

becomes 



Here 



K = -ilim / d^xe'^^0\T(^GG{x)^GG{0)) 10), 
fc-s-o J VSvr Stt / 

^ -z\im f d^xe''^0\T ( ^GG{x)y2zigsqGa^5q{0)] |0). 
Eg. (I273p is obtained by considering 

Lcpv = z^aGG + lidqqGa'y^q 
Stt 2 

and performing path integral. 

Next let us consider duj^^{0) as an 0(0) and use 

d^jf = iq{t^,M}-f,q, 



(273) 

(274) 
(275) 

(276) 



(277) 



where M is the mass matrix of quarks. Then we easily obtain [115 1 

K = m,(0|gg|0), 
K, = ^m, V ^{0\g,qGaq\0). 



2 ^ — ' nig 

q=u,d,s ^ 



(278) 
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So 



oa 



(279) 



leads to vanishing of the s econd term of fl269p 
Finally we obtain 

du + dd 



(0) 
9ttNN 



{N\H^ - Hd\N) 



y-n-NN 



J-IT 



(280) 



Here 



Hq = QsqGaq - mlqq. (281) 

Thus (7^^7v vanishes if we impos e Pe cei-Quinn symmetry. In the absence of Pecei-Quinn 
symmetry, th ere a ppears g^^jq^q 116 |. The contribution of mixing of r] with vr was also 



considered in 



116 |. 



VI. THE EDMS OF MOLECULES 

In this section we consider heteronuclear diatomic molecule which has permanent dipole 
moment. Polar paramagnetic molecules have stronger enhancement factors than param- 
agnetic atoms. Diamagnetic molecules are more sensitive to nuclear P,T violation than 
diamagnetic atoms. 



There are many advantageous points in molecule 117l |. Firstly, the polar molecule is 
polarized by a modest laboratory electric field Eiab but has a vast internal electric field Emt- 
This implies the hugely enhanced stark effect and small fake magnetic field of in 
comparison with atomic case. Secondly, there appears very small energy interval between 
nuclear rotation levels of opposite parity, which is roughly 10~^Ry as will be discussed. Also 
g-factor can be very small etc. 

In general the electric dipole moment D is defined by 



(282) 



where Rj and are coordinates of nucleons and electrons composing molecule. For hetero- 
nuclear molecule 

=< a|D|a (283) 
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and D has the permanent electric dipole moment. 

So the behaviours of heteronuclear molecule and homo nuclear molecule are different. 
First we begin with diatomic molecule with total spin S = case. 
We first give the general rules of diatomic molecule. 

In diatomic molecule, the field has axial symmetry along the two nuclei. Hence the 
projection of L (total orbital angular momentum of electrons) on this axis which is denoted 
by A is conserved. 

The motion of molecule is composed of the orbital motion of electrons, vibrations and 
rotation of nucleus. They interact complicatedly but their interactions are approximated as 
independent motions as the O'th approximation (Born-Oppenheimer approximation) 

^ = ^e^.^r (284) 

and total energy is, therefore, 

E = Ee + E^ + Er. (285) 

They are electronic energy (^ Ry), and vibration and rotation energies of nucleus, respec- 
tively. Let us consider the nuclear motions of diatomic molecule. First we begin with the 
case of total spin (mainly of electrons) S = 0. E^ is considered as a harmonic oscillator and 
its energy are estimated from 

Moola^ ^ . = Ee. (286) 

Here a is the distance between nucleus. M and m are the reduced mass of nuclei and electron 
mass, respectively. Therefore 

Whereas the rotation energy is 

Er = B{K - L)2, (288) 

where K and L are total angular momentum of molecule and electron angular momentum, 
respectively, and 

- 2M^ - 27- (2^9) 
K and the axial component of L are conserved. 

777 

Er = -l{l + l)^-JRy. (290) 
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B 



-J ^ 



•2 



■0 



IITI ]. 



FIG. 19: The vibration {v',v") and rotation (J', J") terms of electron states A and B 



This is much less than the atomic energy interval in general. We show in Fig. [19] the typical 
spectroscopies. 
A-doubling: 



In fl288D K2 and V terms depend on |A|, and K • L oc fi^A ~ (m/Mf^ |ll8 |. 

Hence off diagonal parts are neglected, and +A and —A states are degenerate. 

When we take the relativistic effect into consideration we have another coupling of Spin 
of electrons S (usually nucleon spin can be neglected) with orbital angular momentum of 
electrons Lg and of nucleons L^v. The most important energy shift is A{r)'Le ■ S. 

Selection rule in the electric dipole transition: 

\J' - J\ <1 < J + f (291) 
+ ^-, (292) 



To obtain molecular spectra, we must consider the interactions among the above three terms; 
electron term E^, nuclear vibration i?^, and rotation Er- 

The interaction between Ee and Er is especially important. 
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TABLE II: Angular momenta of diatomic molecule and their projections to molecular axis, n is 
the unit vector of molecular axis. 



angular momentum 


notation 


z component 


electron spin 


S 


S 


electron orbital angular momentum 


L 


A 


nucleus orbital angular momentum 


N 





total angular momentum (without spin) 


K = An + N 


A 


total angular momentum (with spin) 


J = n + N 


n 



First we consider i?e for static nucleus. Unlike atomic case, conserved are not total orbital 
angular moment L and spin S of electrons but their projection to molecular axis 

J, = A + E = n (293) 

which takes the values over A + S,A + S — 1,....,A — S. These states are described as ^^"""^Aq. 
For example, '^Ili/2, "^^3/2 for the states with A = l,S = l/2. 

For atomic fine structure is given by f ll57p . whereas the fine structure for diatomic 
molecule 

dAAT. , . , 
AE = = AA = const. (294) 

We call this relativistic interactions spin- axis interaction, which is composed of spin-orbit, 
spin-spin interactions, as well as the spin and orbital interactions with the rotation of 
molecule. Corresponding to the r elati v e ma gnitudes of these interactions, we can classify 



118| [ll9|. We define the magnitudes of interactions as 



molecule energy levels as follows 
follows. 

LA: the coupling of orbital angular momentum with the axis (the electric interaction 
between the two atoms in the molecule). 

5*^4: the coupling of spin angular momentum with the axis. 
AEf.: the intervals between rotational levels. 

If the distances between terms with different A are larger than both the intervals of fine 
structures (25* -|- 1) and rotational structures, they are further classified into 

• Hund's case a LA > 5^4 > AEr 
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In this case, A, E, Q are well defined and electron state is expressed as "^^^^Aq. For 
the a a transition, 

S' - S = 0, Q' -n = 0, ± 1 (295) 

or 

z/i = AT > z/, = — = AS > z/j = B,{2J + 1) (296) 



Uj{r) = U{r) + A{r)n + B{r){J - L - S)^. (297) 
Here the third term is a perturbation. 

Le and S precess around the internuclear axis z implying that A and S are conserved 
quantum numbers. The total energy is described as 

E = Ee + AeQ + !kj{v + 1/2) + Be{J{J + 1) - 2fi2}. (298) 

In this review we are interested in the transition between parity odd rotation levels of 
the same electron term. 

Hund's case b LA > AEr > SA 

S is not defined. Here the effect of the rotation of the molecule predominates over 
the multiple splitting and total angular momentum J and the sum K = L + N are 
conserved. In this case, S is almost free from molecule (the vector K + S) precessing 
around J, and S is not conserved) 

z/;^ > z/„ (299) 
\K' - K\ <1 < K + K\ (300) 
Ho = He + BK^ (301) 

with K = Az + N and J = K + S. 

UK{r) = U{r) + B{r)K{K + 1) + A(r)A ~ ^ + (302) 

with 

K = A,A + 1,.... (303) 
Here the third term is perturbation. The total energy is 

= f/e + /k^e + BJ<{K + 1) + ^eA^^m^^^l^^. (304) 
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Hund's case c 5*74 > LA > AEr 

Only Vt is well defined. This is the case where the coupling of L with the axis is small 
compared with the spin-orbit coupling. 



Ho = H, + His + 5J2. (305) 

• Hund's case d AEr > LA > S'A 

This is the case where the coupling of L with the axis is small in comparison with the 
intervals in Ej.. 

Ho = H, + BN^ - B{J+r + J-1+). (306) 

• Hund's case e S'A > AEr > LA. 

A. Paramagnetic Molecule 

As we will show, there are a variety of paramagnetic atoms, for instnce, HgF, 
YbF, TIO whose electrons configurations are ™Yb=[Xe]4/i^6s2, ^°B.g=[Xe]4:f^^5d^%s^, 
®^Tl=[Xe]4/^^5(i^°6s^6p^. The selection rules of transitions are 

S' -S = 0, (307) 
A'-A = 0, ± 1 (308) 
E+ ^ E+, S+ ^ S+ for A = 0. (309) 



For BiS molecule 



120|, electron configuration of Bi is [Xe]4/^^5(i^°6s^6p^ and Bi"*""*" has 



one unpaired electron. The electric field of S leads to a mixing of parity odd states: 

\n) = |i/2) = a|si/2, n) + 6bi/2, n) + c|p3/2, (sio) 

Here n = J, = l/2. So 

^ZaYZ\e\de 

^2'2^ = -^ ^(47^-l)a|(iV.W/^ - ^'''^ 

For total J, angular momentum of nuclei takes two values, A'^i = J+1/2 and N2 = J +1/2, 
so the characteristic energy splitting between P-odd states is 

AEr = BN2{N2 + 1) - 5A^i(A^i + 1) = 2B{J + 1/2), (312) 
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which is, for BiS, four to six orders of magnitude smaller than the case of heavy atom. 

2ujdM{^\Hd\(jj) J 



and 



dp 



AEj^, J(J+1) 
(-1)^+1/2^ 



3 X 10^- 



(313) 



(314) 



(J+1/2)(J+1)' 

The effective electric field on the valence electron is proportional to KEmd for polar 
paramagnetic molecule. So it is very advantageous to measure molecular EDM. 



Recently the most stringent upper limit of d^. was reported by using YbF [2J]. Yb 
belongs to the rare-earth elements and its electron configuration is [Xe]-|-4/^^6s^ and Yb"'" 
ion constitutes paramagnetic molecule, f electrons' interaction with the axis of molecule is 
weakened by the deep position of the f electrons and classified as Hund's c class. Their result 
is 



dp 



which sets the upper limit 



:-2.4 ± ^.7 Stat ± 1.5...) X 10-28e cm 



\dp\ < 10.5 X 10"^^ e cm. 



(315) 



(316) 



The other experiment using ThO 



12 1| is also very interesting since a modest laboratory 
electric field Eiat < 100 V/cm fully polarizes a ThO whose internal electric field E^oi is 
lOOGV/cm. (The electron configuration of Thorium is Th=[Rn]6(i^7s^.) This gives another 
advantage for polar molecules. Furthermore, the triplet state ^ Ai of ThO gives the merit of 
g-factor cancellation (see Eq. fl333p ). Also for the other molecules we can expect g- factor 



cancellation, where g- factors are defined by the ratio of spin rotation energy He and iibBz- 
Here He is given by 

He = /3J^ + AS' ■ J - Dn ■ E. (317) 
Here A is the i7-doubling constant and S' is the effective spin and S' = S for Hund's case 



b. The detail of meanings of right-hand side is given in |122| . The expectation value of 
He crosses zero at a specific value of electric field and the molecule becomes insensitive to 
magnetic field at that point. 

One of the problems for molecular EDM is the difficulty of laser cooling compared with 
atomic case. This may be solv ed by first cooling composite atoms and next combining 
them by the Feshbach resonance 123l | and optical trap methods 124l | 125| . The theoretical 
problem is to calculate matrix elements in Dirac-Coulomb + higher order approximation 
(see Appendix ([J])). 
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B. Diamagnetic Molecule 



We will consider TIF as an example of diamagnetic molecule. In searching for molecular 
EDM, we have two tasks. One is to derive dmoi from CP-odd elementary N-N and/or N-e 
interactions. Another is to deduce dp and dn from the observed d^oi- 

The electron configuration of Tl atom is [Xe]4/^^5(i^°6s^6p^ and Tl"^ has a closed electron 
shell. T1+ forms also incomplete shell 6s6p instead of 6s^, 



with 



\n) = |6s, n) + /3\ -^|6pi/2, + y||6p3/2, 



2 Ry a^r(6s, 6p) 
P = —- ^ ^ J = 0.27. 



Eqs — Eq 



(318) 



(319) 



Here Vt = ±1/2, and r(6s,6p) is the radial integral defined by fll80p whose value is 2.3. 
Using dEl]) PI 

Gm^a^ Z^R „ , , , .2 + 7 

p 



{Sl/2\H\pi/2) 



v/27r (iVsiVp)3/2 



where R is the relativistic factor 

R = 



as 



12(27 + 1) V2Zrov4V3 
with To = 1.2fm. As for the nuclear matrix element. 



2-27 



(320) 



(321) 



{k2p X] + hn ^ O-n) ^ k2pj 



(322) 



since a valence proton in Tl atom is si/2- Reference 12|] goes further to get 
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S{Tl) = --ir -ri)d 



(323) 



where r^, are defined by (11170 . From the experimental limit 126 1 



Sexp{Tl) < 0.8 X IQ-^efm^ 



(324) 



we obtain 



dp < lO-^^e cm. 



(325) 
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See (I249p and (12520 for diamagnetic atom. The numerical calculations were estimated along 



the following line of thoughts 
function of TIF is described as 



127j : assuming Born-Oppenheimer approximation, total wave 



*=V'n(rn)V'e(ri)^RK,I). (326) 



Here ipniT^n) describes the motion of Tl nucleus,'?/;e(rj) does F nucleus and electrons with 
respect to the center of mass of Tl nucleus, and iPr^tn, I) the spin and motion of Tl nucleus. 
Let us integrate over ipe and take flB5l) into account. We obtain 



where 

DcT={'ilj^\Y,^n\^n). (328) 

n 

For the present approximation (lB2p . ip^. is given by 

V'e = n,^,(r,) =Tl^Y. ^l<^L{Q^^ (329) 
I 

and so on. 

Anyhow, analytical studies are restricted and we may need more elaborate numerical 
calculations as was done in the case of atomic structures or much more than that case. 
However, it is certain that unknown but very fruitful frontiers are expanding in front of 
this field. Many experiments are preparing or ongoing. In these situations, theoretical 
developments are strongly awaited. 



VII. SUMMARIES AND DISCUSSION 



We have explored the EDMs of quarks, leptons, hadrons, atoms, and molecules. First 
we studied the SM predictions on the EDM and showed that those are far from the present 
experimental upper limits. We have direct signals of new physics beyond the SM from 
neutrino oscillations and muon g-2, and many indirect ones like baryon asymmetry, DM etc. 

Among them, the CP deficiency in baryon asymmetry 77 = ^ ^ 10~^° is especially 
important for searching for new physics. Namely, we can not reproduce r/ via CKM CP 
violating phase only even if we incorporate CP violation due to a ^ term and other radiative 
corrections in the SM framework like GG etc. 



62 



In order to estimate the deviation of phenomena from the SM, we have tried to estimate 
them first in the SM precisely, including the effects of the above mentioned extra terms. 

Next we have explored many theories beyond the SM by focusing on the EDM of elemen- 
tary particles. 

The MSSM and two Higgs doublet model, for instance, give rather large values of EDMs. 
However, those values are mainly due to the ambiguities of the theories themselves. It is 
important to see whether such values are checked to be consistent with the other phenomena 
or not. We think those points are still very insufficient. More predictive models like the 
renormalizable minimal SO(IO) GUT discussed in Section 4.3 often give more stringent 
values which are still several orders smaller than the present upper limits. 

However, the situation is not so pessimistic. Some hope comes from unprecedented col- 
laborations with atomic and molecular physics and elementary particles mainly via brilliant 
developments of laser physics. Most impressive is the new upper limit of the electron EDM 
from polarized molecule YbF. As for paramagnetic atoms, theoretical calculations have been 
developed and seems to be convergent. Whereas, for diamagnetic atoms there are still large 
discrepancies (Table I). Lattice QCD is very promising but it is not convergent in the limit 
of = (Fig. 5). However, it is certain that these situations have been improved rapidly. 
The large parts of such progress have been and will be done by the collaboration of a wide 
field of physics and chemistry. The mutual close relationships among particle, atomic, and 
molecular physics require the wide range of studies over these regions. 

We hope that this review gives some help for these difficult tasks. 

This review is restricted in theoretical part and we have not discussed many excellent 
ideas on the experimental side. The latter is very attractive but is beyond the scope of this 
review simply due to the author's ability. We only briefly explain the mechanism and a list 
of experiments though it is not exhaustive. 

The procedure for the EDM measurement is as follows. First a static external electric 
field E is applied parallel to magnetic field B. The energy splitting is measured as a spin 
precession frequency up^. Next we change E unti-parallel to B whose precession frequency 
is denoted by u^i- Namely, 

/ii/^^ = 2/x-B + 2d-E 

hu-^^ = 2/^-B-2d-E (330) 
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and 

hAu = 4d ■ E. (331) 

Its sensitivity is given by 

, , h 1 1 I , , 
5d = — — 332 

2eKE ^ ' 

Here N: number of sample, r: coherence time, and T: measuring time. K is an enhance factor 

for paramagnetic atoms and molecules given in (11791) and K oc Z^o? . E is a magnitude of 

an internal electric field. So experiments try to get larger values of E, N, r, T. We only 

list up ongoing and planned experiments (see Table 3). We have still more species, solids 

like GGG, Gd2Ga50i2, Gd^Fe^Ou, PbTiO^, Gd^Ga^0i2, solid He, liquid Xe (see Table 

nEDM Collaboration). Please refer to the corresponding sections for the terminologies in the 

experimental features. A few comments are in order. + signature at PbO molecule implies 

the parity under the mirror reflction (reflection under arbitrary plane including molecule 

axis) (see Eq. (l309p ). As for g-factor cancellation in molecular EDMs, ThO and the others' 

cancellation mechanisms are different: the former is due to 

/i = (A + gT.)iiB ^ for A = 2 (333) 

and the latters are due to Eq.f l317p . 

This table is far from being exhaustive but reflects some prospect from a theoretical 
physicist. 

For more detail see, for instance, ECT* Workshop: Violations of Discrete Symmetries in 
Atoms and Nuclei. Nov 15- 19, 2010 |l41 |. 



Finally we will give some comments on the recent results by the Cern Large Hadronic 
Collider (LHC). On Ju l y 20 12 the LHC groups announced the discovery of Higgs-like particle 



around 125 GeV |143| |l44l |. This is not only the discovery of the last unknown particle in 
the SM but also gives the serious impact to the new physics beyond the SM, especially SUSY 
GUT. In this review we have estimated EDM value in the framework of SUSY GUT and 
given large EDM value relative to that of the SM. We briefly explain why 125 GeV Higgs 
mass is serious for SUSY. For tree level Higgs mass satisfies the inequality 

m^<M^|cos(2;g)|, (334) 

with Mz = 91.2GeV, which is obviouly wrong. One loop correction to ruh in CMSSM is 
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TABLE III: A list of ongoing and planned experiments searching for EDM. Superscript * indicates 
estimated sensitivity. 



Species 


Group name 


Features 


muon 






FNAL 


10-^1 e cm* (2015) 




J-PARC 


10"""^^ e cm* (2015) with spin frozen technique 




PSI 


3-4 orders below current limit* (spin frozen technique) 


neutron (all 10^^^ e cm*) 




dn 


ILL (Grenoble) 


\d^\ < 2.9 X 10-2^^ ecm (90% C.L.) [22] 




PSI (Zurich) 


comagnetometer outside of test material 




SNS (Oak Ridge) 


'^He as comagnetometer 




KEK-RCNP (Japan) 


^■^^Xe as comagnetometer [128] 


deuteron 




do 


KVI/BNL/GOSY 


10-^'* e cm* 


paramagnetic Atom 




Cs 


Amherst College 


d(Cs) = (-1.8 =b 6.7 =b 1.8) X 10-^^ ecm [1291 




LBNL 


highly improved magnetic shielding [1301 


Tl 


Berkeley 


< 1-6 X 10-^"(90%Ci..) e cm [20] 


Pr 


GYRIG(Tohoku Univ.) 


K(Pr)=895 EDM measurement starts on 2014 [1311. 


Ra 


KVI (Groningen) 


magneto-optical trap [132] 


diamagnetic atom 




199Hg 


Seattle 


d[^^^^^Hg) < 3.1 X lO-^'* (95% C.L.)[23| 


Ra 


Argonne /KVI 


large enhancement d{Ra) / d{H g) « 10^~^ 


Xe 


@nEDM Collaboration 


polarized liquid Xe droplets 




Tokyo Institute of Technology 
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142| 



where 



M|cos22/3 



3 mf 

4^2 y2 



In^ + t:^ 

2 M| 



= V^h^h^ Xt = At- //cot^, 



12M| 



174GeV 



(335) 



(336) 



with the trihnear Higgs-stop couphng constant At. So 125 GeV indicates heavy stop masses 
and/or large Xt (left-right mixing). Experime ntal serach of SUSY particles also gives large 



sfermion masses, larger than 1 TeV 



145) (l46j. One loop EDM due to the MSSM (Fig.E]) is 



proportional to 0{Mg'^) and heavy Ms reduces EDM. In Fig. [TT] we have assumed Aq = 
since Aq appears at two loop correction in the gauge mediation SUSY breaking and is 
suppressed. In order to preserve rather small Ms and still give large loop correction in 
(13351) . we may take rather large Aq. So we have added the last fourth panel with nonzero Aq 
in Fig. [m More detailed explanations will be given in the review paper of GUT prepared 
by us. 
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Appendix A: SU(6) and Dipole moments 

Both magnetic dipole moment and electric dipole moment are proportional eQcr. 

So we can obtain the information of the ratio of dp/dn = fip/fJ'n from SU(6) in the light 
quark (u,d,s) base if CP violation in the EDM does not affect SU(3) symmetry. They 
are both represented as 

(56|35|56) (Al) 
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Here baryons belong to 56-representation since irreducible representation of qqq = 56 and 
we use that dipole moments are the generator of SU(6). 

\p, 1/2) = ^ {\uud) (2| + +-) - I + -+) - I - ++)) 

+ |mc/m)(2| + -+)-! -++)-! + +-)) (A2) 
+ \duu) (2| - ++) - I + +-) - I + -+))} , 

gaab, 1/2) = ^ l^lnnd) (2| + +-) - | + -+) + | - ++)) 

+^|mmc/)(2| + +-) + ! + -+)-!-++)) (A3) 

-\\uud) (-2| + +-) - I + -+) - I - ++)) 
+cyclic permutations} , 



(p,l/2|ga3b,l/2) 

= 3^ 0(4 + 1 - 1) + ^(4 - 1 + 1) - l(-4 + 1 + 1)) = 1. (A4) 

The corresponding neutron dipole moments are given by the exchange oi u -H- d, resulting 
to I ^ — |. Therefore, 

(n, l/2\Qa3\n, 1/2) 

= 4 + 1 - 1) 4(4 - 1 + 1) + 5(-4 + 1 + 1)) = -| (A5) 

and 

The experimental values of MDM of proton and neutron are [19|] 

/ip = 2.792847356 ± 0.000000023, /i„ = -1.91304273 ± 0.00000045 (A7) 

and the coincidence with SU(6) prediction is good up to quantum corrections. For the EDM, 
compare with the result of lattice calculations FigJSl 
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Appendix B: Multipole expansions 



We will study the multipole expansions of electromagnetic potential = (0, A) due to 
the charged system of finite size. The electric and magnetic fields are defined by 

1 dA 

E = — — grad</), H = rot A. (Bl) 

C CJX) 

Let us assume (as as in the experimental environment) that the electromagnetic field is 
static, that is, the field is time independent. In such case, E (H) is determined only by 
A {(f)). Let us consider a stational motion of chaged particles where Ca charged particles are 
located at and study how the obserber at R feels vector potential A^. 

Here we have neglected the retarded effect of fast particles. If we included it, charged 
distribution has the retarded time dependence and we should replace the arguments as, 

t^t- 1^"^"' , |R-rJ^|R-rJ- '''^^~''"^ (B3) 
c c 

If the scale R r^, (1B2P is expanded around R, 

<^ = ^ - ■ + Y.^ae,ry,dA^ + .... 

= 0(o)+0(i)+0(2) + ____ (B4) 

Then (j)^'"'' is given by 



1 j — z 

^'''^flffrEVsi^eSJ-yje,*). (B5) 



, ■ - + 

m=—l 



where 



= J2^^'a^^'^^rniea,^a). (B6) 

a 

gives electric 2'-pole moment. The superscript (e) indicates electric moment distinguishing 
magnetic counterpart (see (1B25P ). Its continuous representation is 



This comes from 



OO I 1a 

r An 



R- r 



1=0 m=-l 
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First few normalized spherical harmonics Yim are 
Yoo = l/V4^, 

Yw = V3/(47r)cos^, Fiii = vW(8^sin^e^^"^ 



Y20 = V5/(167r)(l -Scos^^), (B9) 



Fa.ii = ±Vl5/(87r)cos^sin^e±*'^, r2,±2 = - Vl5/(327r)sin2^e^2i<^ etc. 
For instance constitute electric dipole moment 

gS=^4, Q'fl, = T^{d.±tdy). (BIO) 
Analogously, vector potential A is expanded as 

Am = I = + 4'^ + A't^ + ... (Bll) 



For instance A) is 



where 



The identity 



leads us to 



4^'^ = (VfeV;-^)r,fcz, (B12) 



Tm = \ I d^rnnMr). (B13) 



j d'rVUr-rrkriJj = (B14) 



where use has been made 

d^mJm = 0. (B16) 

This identity gives 10 constrains. Since Ti^i has 6x3 freedoms and 18 — 10 = 8 physical 
freedoms are remained. We will show that five of eight freedoms constitute M2 moment and 
the remaining three do anapole moment. It goes from subtracting (]B15I) from (1B13P that 



r.. = / d^rr,r„ J. = / d^rr,M„, (BIT) 

where 
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Dividing the Tiki of (1B17P into symmetric and antisymmetric parts w.r.t. l,m, we obtain 



Symmetric part of ( ]i?17l) = -^eifc^M^m (B19) 

o 

with 

Mim = j d^r{riemnr + rmeinr)rnJr- (B20) 
This gives magnetic quadrupole moment. Whereas, 

Anti-symmetric part of (IMTI) = ^ j d^r [5u {rk{rmJ"") - r^Jk) + Ski {J^ - ri(r^J"'))] . 

(B21) 



Here we use the identity obtained from contracting (1B15|) w.r.t. k and / 



j d^r{r^Ji + 2rir^ J™) = (B22) 
Then the anti-symmetric part becomes 

Anti-symmetric part of (\B17h = —{Suak — Skiai), (B23) 

where 

ai = -TT / d^rr^Ji (B24) 



is called anapole moment. 

General expression for magnetic photon corresponding to electric counterpart (]B7P is 



Vim 



1 / An 



-jdh[vx\]-y{r'Yira) (B25) 



/ + 1 V 2/ 

and called 2'-pole magnetic moment (For relativistic case I is replaced by j = |I + s|). 
Appendix C: C,P,T-transformations of Fermi coupling 

We consider the four fermions (current-current) coupling. Here it is concerned with the 
transformation property of Fermion and not with the detailed dynamics, we consider it as 
NONLO'L, where A^, L are spinors and O and O' are a combinations of gamma matrices. 
The most general form are 

GsNNLL + GpN-f^LN-f^L 

+GvNl^NL^''L + GaNi^^I^NL-i^^^L + GrNa^^NLa^" L 
+Gv^lt,NL^^'^^L + Ga^I^.I^NL^^'L (CI) 
+iGs'NNL-f5L + iGp'N-f^NLL + iGT'e^'^'P'^Na^^NLap^L. 
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The first two lines costitute Lorentz scalars, the third line P-odd and the fourth line P,T-odd 
terms. Imaginar i in the last line comes from the Hermiticity of action. The last term of 
the fourth line are also expressed as 

Na^^NLa'^-'^^L or N a ^.^^NLa^'' L (C2) 

since 

e^"'""!, = -il^l^'Yl". (C3) 

C,P,T conjugations are defined by 

C^(t,r) = 7V(t,r), (C4) 
PV^(t,r)=z7Xt,-r), (C5) 
T^\){t, r) = z737V*(-t, r). (C6) 

The fourth line of fICll) is T-odd since 

NN-^-NN, L^^L-^L^^L etc. (C7) 

under T-transformation. 

Appendix D: CP phases in general L-R model and generation number 

Type I (canonical) seesaw is composed of N left-handed and N right-handed neutrino. 
Right-handed neutrino has heavy Majorana mass term. 



LYukau,a = -V^Md^L - -{^^^11^1 + (l^fl)^MRZ/^?} + h.C. (Dl) 

This is described in terms of mass eigen vectors, 



Lm = (iV(i),iV(2)) +h.c. (D2) 

1 yd)* yi2)* j Mn) \ V(^> j ^ ^(2) ' 

2N 2N 

i^iL = Y,UijNjL, uiR = Y,VjiNjR, (D3) 
i=i j=i 

where / = 1, A^. 

So X 2N unitary maqtrices U and V are decomposed into N x N matrices, 
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In the SM 111,^^ = and there is no mixing in neutrino sector. 

For Dirac neutrino case, U^'^^ = V^"^^ = 0. As we mentioned above, there exist (A^ — 
1)(A^ — 2)/2 phases in this case. 

For Majorana neutrino case, U^"^^ = V^^^ = when there exist both left-handed (L-type 
Ni, N]y) and right- handed neutrino (R-type A^at+i, N2n)- In this case 2Nx2N unitary 

V^"^"^ = is added for only left-handed Majorana neutrino case. 



rr _ G 



jlaJl + AjLjI + n (jLj- + ] , (D5) 



where 



jLa = X]/(a;)7a(l -75)t'zL(a;), (D6) 
I 

jRa = ^/(x)7„(l + 75)z/m(x) (D7) 

with / = e,fi,... Mass matrices has rebasing and rephasing symmetries, which does not 
change physics. 

We start with N generation of quarks (Dirac fermions). N x N unitary matrix has 
A^^ real numbers. Of these, 2N-1 is absorved by rephasing of 2N left-handed and right- 
handed quarks. An orthogonal N x N orthogonal matrix has A^(A^ — l)/2 Euler angles. 
The remaining (A^ — 1)(A^ — 2)/2 is the number of phase parameters. Kobayashi-Maskawa 
predicted that there must at least three generations to incorporate CP phase in mass matrix 
2|]. If we relax this arguments to include Majorana neutrino we can use only rephasing of 
charged lepton in MNS mixing matrix whose freedom is N. Therefore, the number of phases 
in MNS matrix is A^^ _ 7v(A^ -i)/2-N = N{N - l)/2. 

If we furthermore generalize the above arguments to include heavy right-handed neutrino 



59|, 



2N 2N 

viL = Y.^i,N,L, yiR = Y,V,iN^R, (D8) 
i=i j=i 

where / = 1, A^. 

So A^ X 2A^ unitary maqtrices U and V are decomposed into N x N matrices, 

U = {U^'\U^'Y , V = {V^'\ V^^^f . (D9) 
In the SM 'my^ = and there is no mixing in neutrino sector. 
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For Dirac neutrino case, f/^^^ = V^'^^ = 0. As we mentioned above, there exist (A^ — 
1)(A^ — 2)/2 phases in this case. 



Appendix E: Expanion in power of 1/c 



Relativistic equation of Fermion in an external electromagnetic foeld is 



{l{p ~ eA) + m) ip = 0. 



(El) 



Let us study the relativistic effect as the deviation from the Schroedinger prescription which 
is obtained by expanding (lEip in power of 1/c. For that purpose we must exclude mc^ from 
energy, which implies to replace ip to ip' 



(E2) 



and 



Substituting 



cct- ( p A ) + f3mc + e$ 



^P'. 



into flE3|) . Dirac spinor is reduced to two component Weyl spinors 



(^^^ - e0 + 2mc^^ x' = c(T ■ (^p- ^A^ x'- 



(E3) 



(E4) 



(E5) 
(E6) 



Retaining only the term 2mc x' in the second equation, we obtain 



= -^(T ■( --A 

2mc V c 



Substituting this into the first equation, we finally obtain the famous Pauli equation. 



(E7) 



The current density is 



1 / 6 \ 2 € 

p A +e<l> — cr-H 

2m V c / 2mc 



j = cilj*aip = c{(j)*ax + X*cr0)- 



(E8) 



(E9) 
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Substituting x of (1E7P into it, we obtain 



2m mc 2m 

In the presence of the EDM (see (l30i) ). j includes pseudo- vector part, 

]d = idN^ X (V'*7V'), 
which in two components approximation is reduced to 

jd = T^V X [0'cr X (p' + p)0], 
2m 

where p and p' are the momenta of and 0', respectively. 
Appendix F: Nonrelativistic approximation 



(ElO) 



(Ell) 



(E12) 



In the heavy nucleon limit, nucleon bilinear forms are approximated as 

N{x)-foN{x) = (5(r), N{x)-fN{x) = 0, 
N{x)-fol5N{x) = 0, N{x)-f-f5N{x) = -cr^(5(r). 



(Fl) 



We are interested in P-odd and T-odd weak interaction in the Fermi coupling between 
electron and nucleons. In the heavy nucleon limit (IFip these interactions are limited in the 
following forms, 

H=-^ [kiNNei-i^e + k2]^e^x^.uNa^xNea^,^ . (F2) 
In the nonrelativistic (heavy nucleon mass) limit it reduces to 



G 

H = i— 5(r)(/i;i7o75 + 4/c2cr ■ 7). 
v2 



(F3) 



Here you should consider H is sandwiched by electron wave functions. For the case of a 
nucleus of charge Z and mass number A, it gives jol 



{Zkip + Nkinhol5 +4:lk2p^(Tp + k2n ^ 



(F4) 



(S1/2I H lpi/2) = g 



Ry 



7 



{Zkip + Nkin) - 8j'^^^{k2p ^ CTp + k2n ^ (Tn) 



8j(l - 7)(A;2p ^ (^np^cTpUp) - ^cTj^ + fcan (^n„(cr„n) - i(T„^ ; 



(F5) 
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The above arguments can be applied for both paramagnetic and diamagnetic atoms. Let us 



apply the above arguments to Cs, Tl, and Xe* atoms 



9j , corresponding to the arguments in 



Section 5.2 in the presence of (lF2p . The wave function for Cs is described as 



|6si/2,F) = |6si/2,F) -3.7 X 10"^^[0.4Uip + 0.59A; 



In 



+ 0.74 X 10"^ 1 



33 



''2p 



\6pi/2,F), 



and, therefore, 



d{Cs) = e{6si/2,F\z\6si/2,F) = -eae x 1.34 x 10 



-10 



0.41A;ip + 0.59fci„ + 0.74 x 10"^ ( F{F + 1) - y ) A;2p 



33 



Here F is the total angular momentum of the atom. The observed value 147l | is 



d{Cs) 



-1.8 ±6.7 ±1.8) X 10"^^ 



e cm. 



For Tl 



For Xe* 



rf(Tl) = ettB ■ 0.96 X 10-^(0.4fcip ± 0.6fci„ - 2 ■ lO'^Ajap) 
d{Xe*) = -1.3 ■ 10~i°eaij(0.41A;ip ± 0.59A:i„). 



Appendix G: Strong CP violation 



In the QCD world, the true vacuum is described by the 9 vacuum, 



\0) = Y^ e"™^ \n), {n = integer) . 



|<i^|„_„,e'/"'" 



-iud 



[dAU 



i f Ldf^x 



Using that An gives 



n 



(F6) 



(F7) 



(F8) 

(F9) 
(FIO) 



(Gl) 



(G2) 



(G3) 
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and substituting (1G3P into (lG2p we obtain 



with 



y\e-'"'\e) = J2j [d^W^^^" (G4) 



Le// = L + j^ I d'^xTr ( G^.G'^'^ ) (G5) 



Appendix H: U(l) problem 

6 term in (120 ip comes from the fact that the vacuum in QCD is 6), whereas GG term in 
(12161) does from quantum anomaly, occurring irrelevant to Abelian and non Abelian. In this 
appendix we will show that these two terms are closely related and lead us to solve U(l) 
problem. 

In the following discussions we consider mass zero quark limit, and Nj = 3, up, down, 
strange quarks. Chiral invariant action has originally Ul{3) ® f/i?(3) symmetry. If, as we 
have considered, QCD vacuum is quark condensate 

{uu) = (dd) = (ss), (HI) 

action symmetry is reduced to the flavor symmetry U(3) and generates 3^ NG bosons. They 
are ti"^, 7r°, K"^, K^, , r]s, and Here the first eight particles constitute octet and the 
last a singlet. The observed mass eigen states, rj and t]' particles, are the linear combina tions 



of rjs and ?7o, and their masses are = 550MeV, m^/ = 958 MeV. Weinberg showed 148 1 
that the observed m^/ is too heavy for predicted NG boson, 

m^/ < y/SniT,. (H2) 

This is one of the U(l) problems. Ano ther is concerned with rj tt+tt^tt^ process. 
Let us explain these problems 



149] : The octet axial vector currents satisfy 
d^JS^ = famlr (a = l,...,8) (H3) 

and 

Sai,mlf! = I'^^^^zK [ d'xe-^'%Q\T{d^Jt,m^Jl{x)) |0) 

iK I rf^xe-^^-(0|T {d^J^^iO)Jl{x)) |0) (H4) 



.mi — 



ml 



+ I d'xe~^'^0\6{xo) [d'-J^.iO), JUx)] |0) 
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In the low energy limit, if there is no massless pole, this reduces to 



W/a = ^ J d'x{0\6{x'M,{x),d^J^^{x)]\0) (H5) 

Whereas, isosinglet axial vector current constitute ABJ anomaly (I216p . The isosinglet 
can be described as a sum of SU(3) octet and singlet, 



■^sm — ^"^5 J + \/ ^Jbix 



(0) 



(H6) 



with 



4? = -^{u-ff,-f5U + mdd-ff,-f5d-2s-ff,-f5s). 



J, 



(0) 
5fi 



(H7) 
(H8) 



Taking fl216p . f l218p . and f l219p into considerations, we obtain the same equation for isosinglet 
case as f lHSp by replacing J^^ with Js^, 



2 £2 



.ml — k'^ 



rrin 



iK I d^xe-'^'' {Q\T [d^" Jl{Q)Jl{x)] |0) 



d^xe-'^''{Q\5{xQ) d^^rJO), |0) 



(H9) 



Here /o is the isoscalar meson decay constant. If there is no zero mass pole like the octet 
cases, this relation is same as the octet case flHSp except for J^^ replaced by Js^ , and we 
obtain 

mlfo = <fl (HIO) 
So, if SU(3) is good symmetry, it goes from ( 1H6P and ( IHlOp that 



/o > ^/., 



(Hll) 



ess particle couples to J^, then first term 



which directly leads to (IH2p . However if any mass' 
of ( IHSp does not vanish and we can evade ( IHIOI) 150| . 't H ooft showed that this is indeed 



the case if we take 9 vacuum into consideration correctly 
solution compatible with quark condensate |l52| : 



15 1| . Also Witten proposed a 



(H12) 
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where 

'd^En\ 1 / 1 



j d^x{T {Tr{G{x)G{x))Tr{G{Q)G{Q))^) (H13) 



(vr+TT-vr^lr/) = ^" lim(7r+7r-|g^J5^(fc)|r/). (H14) 

I'^TTlq k^o 

The right-hand side vanishes due to momentum conservation. However, it is experimentally 
obse rved as r(?7 — )■ vr+vr"??^) ^ 200eV. This process is occurred via SU(2) violating operator 



153 | 



and 



where 



C = -{rriu — rridjiuu — dd) (H15) 



(37r|£|7/) ^ -j^^ , (H16) 



A = (7r7r|(m„u75U + md(i75(i)|?7) 



1 



^.{nn\d'^J',iOM- (HIT) 

So this process is suppressed by axial vector current conservation even if mu ^ md- This is 
another U(l) problem. 



Appendix I: Schiff moment 



There are several origins for the Schiff moment. Here we discuss the Schiff moment 
induced by the nuclear EDM when the charge and the EDM distributions (pg and p^, re- 
spectively) in the nucleus are different [l2|. 

The interaction of the electron with the dipole moment of finite size nucleus is 



Vs = j dV[p,(r')-p,(r']d^- V- 



— e 



r — r' 



K = j d'r'[pd{v') - PqiAdN/^n^iV^vJ^. 

Here we may assume 154l |: 
Pq is spherically symmetric. 

d]\f coincides with the EDM of valence nucleon, = dp^n<^. 
Pd is due to the valence nucleon. 



(II) 
(12) 
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Then 



mn 



1 



where n = r'/r'. Let us divide V/VmVn as 



5 



(13) 



(14) 



The first term corresponds to the electron interaction with the 2^-pole moment of the nucleus. 



The second term 



mn 



m I 

5 



■ V 



Here we use 



154| 



(15) 



/Stt 



cr ■ nn = -a 

3 



Att[Yo (g) a]^i 



(1) 



2cr ■ nn + cr = V47r ( - [Fq ® cr]^^^ 



2^2 



(1) 



(16) 



where 



•nil ,ms 



with the Clebsch-Gordan coefficient {l,mi; ^,ms\j,m) related with 3j-symbol 



(17) 



{ki,qi,k2,q2\K,Q) 



V2K+1 



ki k2 K 
qi q2 -Q 



(18) 



The following equation is the Wigner-Eckart theorem (the definition of reduced matrix 
element) ( || || ), 

1 



{Km\0 jM\i^fn') 



V2J + T 



{jmJM\jm'){K\\Oj\\K) 



(19) 



The following arguments are indebted to discussions with T. Sato 
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where k is defined by (11831) and 



|k) = \[Yi{n)^x]ij))- 



(110) 



It should be noted that the reduced matrix has no dependence on m, m'. M. 
For J = 1 case 



{k\0\k) = {k\J\k) 



{k\\0\\k) 



{k\\J\\k) 
(fi:||J||ft:) = VJ(J+1)(2J+1) 



,^\\k) = 2|«:|(-l)(-l)l'^l(j,i;j,-i|l,0) 
1 - 2k 



V3 
2(1 + 

1 / 2j + l 



2'-^' 2' 
for / = 



for / = 2 



2V J(J + 1) 



X 



1 - 2fi: for / = 
-V2{1 + k) for/ = 2 



Thus we obtain 



-Pain ) pJk, ' 



(j)-v 



So 



r^-(/s; ) — r^-(K — -) 



jV 

jU + 1) 



47r5(r) 



where the mean squared radii was defined by 



(111) 
(112) 



(113) 



(114) 



(115) 



(116) 



(117) 



In this derivation we assumed that the nuclear charge is uniformly distributed over a sphere 
get the final expression for the Schiff moment. 



of radius rn = 1.2 x 10 ^^A^^^ cm, and r1 = |rn. Also we may assume r'^ = r'^ Q]. Then we 



''-^^•'^"°25j(j + l)" 



(118) 
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Appendix J: Effective Hamiltonian in molecule 

We have said that there appears huge internal electric field Ej„f in polar molecule. Here 
we consider how to estimate -Emi- The Dirac-Coulomb Hamiltonian is 

Ho = ^{ctti ■ Pi + f3imc'^ + Vnuci{ri)} + — (Jl) 
and P,T-odd perturbation (the intrinsic part of HpTv) is 

ij' = -4$^Ao--Er*. (j2) 

i 

with 

Ei,int = - V J Vnuci (r^) + 2 XI — ) • ( J3) 



Here electric field is given by Eq. flBip with = X]j{^ri«d(rj)}+Xli<,- ir-}- We are considering 



c dt 



a static field, = 0, and H' is represented as 



ff' = 45^[A'T.-V„ i/-T] (J4) 

i 

Here T is kinetic term of electron 

T = J]{cai-p, + A^c2}- (J5) 

i 

The expectation value w.r.t. the eigen function of gives 

V„ i7]|*) = (J6) 

Whereas, 

X [AcTi ■ Vi, T] = i X ^{[/JiCri ■ Pi, q;j ■ pj] + [/JjCTj ■ p^, - l)mjc]} 



Here the second term vanishes and the first term gives 

(X 2/3i75p- for % = j. 
for z 7^ j. 



So 



^=1^1 (J7) 
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Thus H' of flJ4l) is rewritten as 



155 | 



i 

and we obtain finally 

2zc(V^o|/375pVo) = -4cp2!3((^tx)- 
Thee enhancement factor is given by 



E — E„ 



+ h.c. 



(J8) 



(J9) 



(JIO) 



So the detailed calculations are reduced to the electron wave functions in atoms and 
molecules. For molecular case, unfortunately, only H2 can be solved in the Born- 
Oppenheimer approximation 156| . However, its perturvation expansion around at omic level 
is also interesting since this method is applicable to t he other diatomic molecule 157l |. For 
more detailed explanation for diatomic case, see 1581]. 
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